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Division, Navigation and Guidance Research Branch, the NASA Langley Research 
Center (LaRC), Hampton, Virginia. Mr. Henry J. E. Reid, Jr. served as Technical 
Monitor on the contract. 

This is a final technical report which documents the results of research per- 
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Modification 5 to the original contract. 

The effort was directed by Mr. John Zvara, President and Technical Director 
of ASIo Mr. William C. Hoffman served as Project Engineer. Dr. Walter M. Hollister 
of the Department of Aeronautics and Astronautics at the Massachusetts Institute of 
Technology (MIT) contributed to the study as technical consultant and co-investigator. 
Dr. Robert W. Simpson, Director of the MIT Flight Transportation Laboratory and Dr. 
Arthur E. Bryson, Jr., Chairman of the Department of Aeronautics and Astronautics, 
Stanford University, also served as technical consultants. 


lii 


AEROSPACE SYSTEMS, INC. • ONE VINE BROOK PARK • BURLINGTON, MASBAOHUSETTS 01B03 • (BIT) 376-7017 


TABLE OF CONTENTS 


Section Page 

FOREWORD ■ iii 

LIST OF ILLUSTRATIONS jx 

LIST OF TABLES xiil 

1 INTRODUCTION ^ 

2 OPERATIONAL CONSIDERATIONS 5 

2.1 VTOL All-Weather Operations 5 

2.2 Basic Spiral Parameters 8 

3 SPIRAL GUIDANCE FORMULATION 15 

3 . 1 Basic Concept 15 

3.2 Equations of Motion 15 

3.3 Spiral Angle as Independent Variable 17 

3.4 Time in Spiral for Arbitrary Turn Angle 18 

3.5 Nominal Spiral Specification 72 

3.6 Wind Effects on Nominal Spiral 73 

3.7 Feedback Guidance Laws 25 

4 PERTURBATION GUIDANCE ANALYSIS 27 

4.1 Simplified 3D Guidance Without Wind 27 

4.1.1 Horizontal Plane Guidance 28 

4.1.2 Vertical Plane Guidance 30 

4.1.3 Example 31 

PRECEDING PAGE BLANK NOT FILMEQI 

V/ 


aerospace systems, INC. • ONE VINE BROOK PARK . BURLINGTON, MASSACHUSETTS OIBOB • (S'!?) S7S-7B17 


TABLE OF CONTENTS (Continued) 


Section Page 

4.2 4-D Guidance 32 

4.3 Guidance Gains 42 

4.4 Error Sensitivity 45 

5 WIND ESTIMATION 47 

5.1 Kalman Wind Estimator Formulation 47 

5.2 Simplified Wind Estimator 50 

5.3 Wind Estimator Dynamics 55 

6 PHASE PLANE ANALYSIS 59 

6.1 The Phase Plane 59 

6.2 Guaranteed Spiral Capture 51 

6.2.1 Right Hand Turns 51 

6.2.2 Left Hand Turns 64 

7 SPIRAL TRANSITION SEGMENTS 57 

7.1 Nominal Path for Entry Transition to Spiral 57 

7.1.1 Right-Hand Spiral Entry * 57 

7.1.2 Left-Hand Spiral Entry 59 

7.2 Nominal Path for Exit Transition from Spiral 70 

7.2.1 Right-Hand Spiral Exit 70 

7.2.2 Left-Hand Spiral Exit 72 

7.3 Transition Segment Guidance Law 72 

7.4 Hover and Touchdown Phase 74 

vi 


AEROSP ACE 8YBTE!VI8| INC. • one VINE brook pabK • Burlington, masbachubettb 0 iB 03 • (ai 7 ) 27B-7B17 


TABLE OF CONTENTS (Continued) 


Section Page 

8 SIMULATION., 77 

8.1 Interactive Simulation 77 

8.2 Program VALT Simulation 78 

8.2.1 Initialization Calculations 80 

8.2.2 Calculating Nominal Parameters 83 

8.2.3 Wind Estimates 88 

8.2.4 Guidance Laws 89 

8.2.5 Measurements 90 

8.3 Spiral Approach Example 92 

9 SUMMARY AND CONCLUSIONS 109 

10 RECOMMENDATIONS FOR ADDITIONAL RESEARCH 113 

10.1 Spiral Guidance Flight Test Program 113 

10.1.1 Purpose of Test 113 

10.1.2 Test Equipment and Facilities 113 

10.1.3 Flight Procedures 114 

10.1.4 Data Obtained 114 

10.1.5 Data Analysis 115 

10.2 Guidance for Spiral Exit and Landing 115 

10.3 Wind Estimation in the VALT Flight program 116 

REFERENCES 119 

vM 


ASROSPACE BYBTEIVIEi INC. • one Vine brook park • BURL.iNaTON, massachusette oibos • (ai7) &7e-7Si7 


LIST OF ILLUSTRATIONS 


Figure No. 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 


Noninterfering VTOL Approach to Conventional Airport.... 

Example Spiral Descent for JFK Inf I Airport 

Spiral Descent Approach to CTOL Airport 

Spiral Descent Approach to Heliport 

Nominal Spiral Parameters - No Wind 

Effect of Wind on Spiral with Constant Airspeed 

Block Diagram of Spiral Guidance Concept 

Nomenclature for Spiral Descent .‘ 

Nomenclature for Determining t* (0 ) 

Nomenclature for Spiral Descent Guidance Scheme 

Linearized Block Diagram of Horizontal Guidance System .. 

Spiral Guidance Response for Various Sets of Gains 

Effect of Finite Bank Rate 

Nomenclature for Wind Analysis 

Estimated Wind Speed with Simplified Estimator 

Estimated Wind Direction with Simplified Estimator 

Radial Position Accuracy with Simplified Wind Estimator ••• 

Linear Signal Flow Diagram for Wind Estimator 

Wind Estimator in a Steady Shear 

Linear Signal Flow Diagram for Wind Estimator During 
T ransition Segment 

Phase-Plane Plot of Guidance Law * 

Phase-Plane Performance of Horizontal Spiral Guidance Law 


Page 

6 

7 

9 

10 

11 

12 

15 

16 
19 
27 
33 

43 

44 
47 

52 

53 

54 

56 

57 

58 

59 

60 


preceding page blank not EHJlEn 


AEROSPACE SYSTEIVIBi |NC« • one vine brook park • BURLiNaroN, masbachusettb aieo3 • (ei7) S7 s-7B17 


LIST OF ILLUSTRATIONS (Continued) 


Figure No. Page 

23 Phase-Plane Examples 62 

24 Modified Guidance Law 63 

25 Example of Modified Guidance Law for = 60° and 90°. . . 64 

26 Independent Variables for Straight Segments and Circular Arcs. . . 67 

27 Entry Transition to Right-Hand Spiral 68 

28 Entry Transition to Left-Hand Spiral 69 

29 Exit Transition from Spiral 70 

30 Modified Transition Segment Guidance in Phase Plane 73 

31 Pad-Centered Frame for Hover and Touchdown 74 

32 Overall Flow Diagram of Program VALT 79 

33 Segments of Nominal Spiral 80 

34 Overall Flow Diagram of Spiral Guidance Algorithm 31 

35 Rotated Coordinate Frame for Spiral Entry 84 

36 P^otated Coordinate Frame for Spiral Exit 36 

37 Example Spiral Approaches to JFK 93 

38 Nomina! Spiral Approach to JFK - Position Error Histories 95 

39 Nominal Spiral Approach to JFK - Velocity Error Histories ...... 96 

40 Initial Crosstrack Error of -1 nm 99 

41 Initial Along-Track Error of 1 nm 100 

42 Initial Groundspeed Error of 20 kts 101 

43 Initial Altitude Error of -300 ft 102 

44 Initial Ground Track Error of -30 Deg 103 

45 Constant Wind Direction Error of -90 Deg 104 

X 


AlSHOSPAOE 8Y8TEMG1 INC. • one VINE SROdK PARK • BURUNBTON, MASBACHUBETTB OIBOa • (b 17 ) 


LIST OF ILLUSTRATIONS (ConMnued) 


Figure No, Page 

46 Wind DirecHon Shear Error of -90 Deg/100 Ft 105 

47 Nominal Wind Speed Error of -15 kf 106 

48 Nominal Wind Direction Error of -90 Deg 107 

49 Possible Heliport Locations Relative to the Descent Spiral 115 


xi 


aerospace 8YBTEM8, IIMC. • one vine brook park • Burlington, masbacHUsettb oieo3 • 


(B17) 878-7817 



LIST OF TABLES 


Table No. Page 

1 Cons^ant• Groundspeed Versus Constant Airspeed 13 

2 Operational Procedures for 3-D and 4-D Guidance 14 

3 Approximate Solution to Elliptic Integral of the Second Kind 21 

4 Wind Effects on 360° Turns (V = 60 kt, r = 2000 ft) 24 

5 Nominal Turn With Bank Control Only 35 

6 Effect of Initial Radial Error With Bank Control Only 35 

7 Effects of Initial Heading Error With Bank Control Only 36 

8 Effects of Airspeed Error With Bank Control Only 36 

9 Effect of 5-Knot Wind With Bank Control Only (90°) 37 

10 Effect of 5- Knot Wind With Bank Control Only (180°) 37 

11 Effect of 20-Knot Wind With Bank Control Only 38 

12 Effect of Including Wind Estimate in Nominal Calculation 38 

13 Effects of Airspeed Error With Coupled Bank & Throttle Control .... 40 

14 Effects of Airspeed Error With Uncoupled Bank & Throttle 

Control 40 

15 Effects of Wind With Coupled Bank & Throttle Control 41 

16 Effects of Wind With Uncoupled Bank & Throttle Control 41 

17 Recommended Values for Guidance Gains 91 

18 Spiral Guidance Parametric Analyses 97 

19 Results of Parametric Analyses 98 

PRECEDING PAGE PLANE NOT gTLMlSa 


xiii 


AEROBPACE BYETEMS, INC. • one vine bhook park • burungton, Massachusetts 01B03 . (ai 7 ) 272-7017 


SECTION 1 
INTRODUCTION 

VerHcal fakeoff and landing (VTOL) systems have been recommended by 
many experts in the aeronautical and transportation fields, including members of the 
President's Scientific Advisory Committee, as the logical form of transportation for 
the high-density, short-haul transportation markets in the 1980s. Such systems span 
a broad range of operations, from the intra-urban to the inter-city shuttles operating 
between downtown or nearby vertiports as well as to conventional airports. 

Before a viable VTOL system can become a reality, technology developments 
are needed in several areas. For example, the technology for an economical, 150- 
passenger class VTOL with reasonably high cruise speed and acceptable passenger ride 
qualities for the inter-city market must be developed. At the other end of the spectrum, - 
advanced helicopter development is needed to improve ride comfort and reduce main- 
tenance requirements for the very short-haul, medium-density market. These vehicle 
design areas are receiving considerable attention in various programs sponsored by 
NASA, the U.S. Army, and the aircraft industry. However, to effectively utilize 
these vehicles, 'ind to exploit their unique characteristics for minimizing noise and 
both air and ground space requirements, corresponding odvances must be made in 
handling qualities, operating procedures and techniques, and avionics. 

The NASA Langley Research Center (LaRC) has undertaken a research program 
to develop the navigation, guidance, control, and flight management technology base 
needed by Government and industry in establishing systems design concepts and operat- 
ing procedures for VTOL short-haul transportation systems in the 1980s time period. 

The VALT (VTOL Automatic Landing Technology) program encompasses the investigation 
of operating systems and piloting techniques associated with VTOL operations under all- 
weather conditions from downtown vertiports; the definition of terminal air traffic and 
airspace requirements; and the development of avionics including navigation, guidance, 
controls, and displays for automated takeoff, cruise, and landing operations. 

Aerospace Systems, Inc. (ASI) is conducting a research effort for LaRC in 
support of the VALT program. In a previous report (Reference 1), ASI analyzed the 
navigation and guidance requirements for commercial VTOL operations in the takeoff, 
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cruise, terminal area, and landing phases of flight in weather conditions up to and 
including Category III. The study was limited to two types of rotorcraft vehicles 
(pure helicopter and compound helicopter) and three types of services (intra-urban, 
inter-city, and conventional airports) were considered. Applicable navigation tech- 
nology and systems (such as Omega, Loran, inertial, and microwave landing systems) 
were examined to define present system shortcomings, to Identify areas where technology 
advances are required, and to select candidate systems and conceptual approaches. A 
multi-configuration "straw-man" system design was prepared, and representative opera- 
tional procedures and trajectories defined. A limited flight eva fixation program was 
conducted to investigate VTOL operational procedures and current ncsvigation systems 
and to verify analytical results. A comprehensive digital computer simulation (Program 
VALT) was developed to provide a means for evaluation of VTOL guidance and navi- 
gation system performance. 

One of the conclusions reached in Reference 1 was that curved decelerating 
approaches will be required for safe, efficient, and independent VTOL operations. To 
facilitate these maneuvers, a spiral descent technique was proposed as a possible stan- 
dard VTOL approach procedure. The spiral descent uses minimal airspace, accommodates 
arrivals from any direction, and can service multipad landings. The spiral approach also 
provides the benefits of a vertical descent, but avoids the vortex ring state, maintains 
a stable airspeed, and uses less fuel . Flight evaluations were conducted for several 
spiral descents, and demonstrated the concept feasibility in terms of ride quality, 
vehicle capability and pilot workload. 

As a result. Reference 1 recommended that additional research should be con- 
ducted to define the navigation and guidance system requirements for IFR VTOL spiral 
descents in the presence \of winds. It suggested a study effort to establish recommended 
values for airspeed, bank angle, descent rate and protected airspace; to formulate a 
number of feasible guidance laws for spiraling flight; ond to develop filtering algorithms 
for navigation during this phase. The guidance and navigation algorithms could be 
evaluated with the simulation Program VALT. The necessary commands to display on a 
flight director and/or the signals to feed an automatic pilot should be established. 
Finally, one or more spiral guidance laws should be recommended for further evaluation 
in the LaRC VTOL real-time simulation, and eventual implementation in the VTOL 
flight test program. 
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Subsequently, NASA LaRC decided to pursue the spiral guidance approach 
concept with additional research. The present report describes the results of ASI‘s 
investigation of the guidance and navigation req^iv^ments for VTOL spiral descents 
in the presence of winds. Models have been developed to describe the spiral maneuver 
and candidate guidance laws have been formulated and analyzed. An important ele- 
ment of the guidance scheme is a unique wind estimator which uses the perturbations 
in bank angle and heading to improve the knowledge of the winds. Finally, recommen- 
dations for additional research, including a flight program, have been outlined to 
evaluate the spiral guidance concept. 
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SECTION 2 

OPERATIONAL CONSIDERATIONS 

2.1 VTOL ALL-WEATHER OPERATIONS 

The achievement of efficient all-weather VTOL commercial service to con- 
ventional airports requires that VTOL aircraft be operated essentially independent of 
existing CTOL aircraft route structures and procedures. The establishment of such 
independent VTOL operations during the takeoff and landing phases must take into 
account noise restrictions^ obstacle clearance, VTOL aircraft performance capabilities, 
and Interaction with CTOL air traffic. Previous studies of navigation and guidance 
requirements for VIOL operations at CTOL airports (Reference 1) have suggested the 
use of spiral descents by VTOL aircraft as a means of achieving efficient, safe opera- 
tions that are independent of CTOL aircraft route structures and procedures. 

At the present time, IFR operations of CTOL aircraft create a wall of airspace 
along the active runway that is typically 10 miles long, 1500 feet high and some hun- 
dreds of feet wide. The protected airspace is utilized by CTOL aircraft making ILS 
approaches, departures and missed approaches; it creates a problem for VTOL traffic 
desiring to cross the active runway with no interference to the CTOL traffic flow. A 
suggested noninterfering VTOL approach to a conventional airport is illustrated in 
Figure 1 . Under this concept the airspace directly over the active runway is partitioned 
such that the CTOL traffic remains below 500 feet AGL over the runway, and VTOL 
traffic has free access to cross over the runway perpendicular to the CTOL traffic at 
1000 feet AGL. Failure to allow IFR crossing of the CTOL runway in this way would 
necessitate on approach capability to both sides of each CTOL IFR approach course, 
followed by an air taxi across active runways. Such operations across active runways 
would have to be conducted under CTOL air traffic control, which could cause exces- 
sive delays and would make the VTOL traffic dependent on the CTOL operations. 

After the 1000 feet crossing, the VTOL will take 2 minutes to descend at a 
rate of 500 ft/min. To expedite the approach, it is desirable to begin the turn to the 
pad as soon as possible. Using a standard turn rate (3 deg/sec) at an airspeed of 60 
knots, the approach will consist of a spiral descent with turn radius of about 2000 feet, 
which fits conveniently into the typically available airspace. For example. Figure 2 
illustrates a 2000 ft spiral superimposed on a diagram of New York's Kennedy 
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VTOL APPROACH 



Figure 1 Noninterfering VTOL Approach to Conventional Airport. 
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Figure 2. Example Spiral Descent for Kennedy International Airport (JFK). 



InternaHonal Airport. This solution has the advantage of keeping the two traffic flows 
independent and without mutual interference. It also influences the characteristics of 
the VTOL approach procedure and the associated guidance requirements. 

This spiral technique can be generalized to accommodate arrivals from all 
directions and could handle multiple helicopters in a "descent tube," as shown in 
Figure 3 . The descent tube would be established in a sector of the airspace not used 
by CTOL aircraft, with rotorcraft entries occurring above the CTOL approach/departure 
patterns. The spiral is not restricted to CTOL airports, but can be applied to single or 
multiple pad V-ports as well, with provisions for a missed approach if required 
(Figure 4). . 

The spiral descent retains most of the advantages of a vertical descent, but 
requires less power, maintains forward airspeed and controllability, and avoids the 
vortex ring state. The problem with the spiral descent operation is that it presents a 
difficult guidance and navigation task, particularly in a wind. The steady-state turn 
introduces an additional integration into the control loop, making it more difficult to 
stabilize. The maneuvering aircraft may have difficulty in obtaining accurate naviga- 
tion information at a sufficient data rate. The present study was undertaken to evaluate 
such problems and to formulate guidance concepts for spiral descents during the planned 
NASA VALT flight test program. 

2.2 BASIC SPIRAL PARAMETERS 

The design of the spiral descent in the vertical plane is quite simple, requiring 
only an appropriate selection of the vertical speed. However, in the horizontal plane, a 
number of parameters must be properly chosen to satisfy airspace limitations, aircraft 
constraints, passenger comfort, etc. Figure 5 illustrates the principal spiral param- 
eters in the horizontal plane for the case of no wind. Without wind, the airspeed V 
is equal to the groundspeed, and the heading angle ^ is perpendicular to the spiral 
angle 9. The heading rate is 

S-. tan 0 (1) 

V 

where g is the acceleration of gravity. 

The bank angle 0 required to maintain the steady spiral descent is given by 
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Figure 3. Spiral Descent Approach to CTOL Airport. 
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NORTH 



Figure 5. Nomina! Spiral Parameters - No Wind. 

gtan0 = vVr (2) 

Thus, from Equations (1) and (2), specifying any two of the four spiral parameters 
determines the remaining two. For example, a 60 kt spiral with a 2000 ft radius requires 
a bank angle of 9.05° and a heading rate of 2.9°/sec. 

The effect of a wind on the spiral descent is to produce a continuous variation 

in the spiral parameters. The airspeed V and groundspeed V now differ, which also 

results in changing heading rate, bank angle and crab angle (between V and V ). This 

9 

is depicted in Figure 6 for a spiral descent which maintains constant airspeed in a 
steady wind. The magnitude of the crab angle is greatest with a direct crosswind 
(points 1 and 3), while the bank angle is about the same as the zero-wind case. The 
highest groundspeed and maximum bank angle, occur downwind (point 2), with the 
reverse situation at the upwind position (point 4); the crab angle is zero when the 
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WIND 



Figure 6, Effect of Wind on Spiral with Constant Airspeed. 
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aircraft is directly downwind or upwind. If a constant rate of descent is maintained, the 
changing groundspeed will produce a nonlinear variation in altitude with spiral angle. 

The preceding discussion has implied two philosophies in flying the spiral 
descent in a wind: constant airspeed or constant groundspeed. In the previous analyses 
of commercial VTOL operations, the groundspeed has been specified to simplify 4- 
dimensional navigation. However, several factors favor constant airspeed for the 
spiral approach: the descent airspeed can be chosen for fuel economy; large airspeed 
variations may take the rotorcraft Into unsafe operating regions; and continually- 
varying longitudinal accelerations would be required to maintain constant groundspeed 
while spiraling in a steady wind. Some consequences of constant groundspeed versus 
constant airspeed in a steady wind are summarized in Table 1 . 

Table 1. Constant Groundspeed Versus Constant Airspeed . 



CONSTANT 

GROUNDSPEED 

CONSTANT 

AIRSPEED 

Acceleration Vector 

Constant magnitude 
(directed to center of spiral) 

Variable magnitude 
(directed perpendicular 
to airspeed vector) 

Longitudinal Acceleration 

Variable 

Zero 

Airspeed 

Variable 

Constant 

Heading Rate 

Variable 

Variable 


Unfortunately, the constant airspeed approach complicates the timing problem 
of precise inter-aircraft spacing for 4-D guidance. For 4-D guidance at constant air- 
speed, the wind must be accurately estimated to predict the groundspeed. However, 
the wind knowledge is also very important for 3-D guidance as well, and a wind 

estimator will be required in any event. A comparison of important operational pro- 
cedures for 3-D and 4-D spiral approaches is shown in Table 2 . 
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Table 2 , Operafional Procedures for 3-D and 4-D Guidance. 


Approach 

Plate 

Information 

Definition of specific spiral approach procedure 

• Center and radius of spiral tube 

• Entry and departure angles 

• Altitude versus angle 

• Zero-wind time 

Approach 

Clearance 

3-D 

ATC provides wind estimate 
Clearance assumes constant 
airspeed during spiral 
Clearance assumes constant rate 
of descent b'jsed on nominal wind 

4-D 

ATC provides wind estimate 
ATC directs time to exit spiral 
based on constant airspeed and 
wind estimate 
Clearance assumes altitude 
versus time profile based on 
constant rate of descent and 
nominal wind 

Descent 

Procedures 

3-D 

VTOL maintains radius with bank 
angle 

Actual wind is estimated to improve 
nominal control 

Airspeed is maintained constant 

Descent rate is adjusted to achieve 
final altitude at final angle 

4-D 

VTOL maintains radius with 
bank angle 

Actual wind is estimated to 
improve nominal control 
Airspeed is adjusted to achieve 
nominal exit time 
Descent rate is adjusted to 
follow altitude versus time 
profile 
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SECTION 3 

SPIRAL GUIDANCE FORMULATION 
3.1 BASIC CONCEPT 

The basic concept followed in the spiral descent scheme is that of guiding 
the rotorcraft to a nominal path which can be prespecified for any specific landing site. 

The genera! concept is illustrated in Figure 7 i n block diagram form. The navigation 
system takes measurements of the actual rotorcraft state x, and produces estimates of 

A A 

the state vector x, and of the independent variable t . The latter is used to calculate nominal 
values of both the state x* and the control u*. The nominal and estimated states are 
compared, and any errors are fed back through the guidance gain matrix to generate 
control corrections Au. These corrections are subtracted from the nominal control to 
obtain the actual control u. Uncertainties enter the system as wind disturbances on 
the rotorcraft, navigation system inaccuracies, and approximations used in the nominal 
calculations in the guidance gains. 



Figure 7. Block Diagram of Spiral Guidance Concept. 
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3.2 


EQUATIONS OF MOTION 

The basic dynamics of the spiral descent involve five state variables (Figure 8): 


r = rotorcraft radial distance from nominal spiral center 
0 = rotorcraft azimuth around spiral, clockwise from North 
^ = heading angle ofoirspeed vector, clockwise from North 
V = rotorcraft airspeed 
h = rotorcraft altitude 


NORTH 



Figure 8. Nomenclature for Spiral Descent. 

« 
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The basic differential equations describing the actual spiral descent are: 


V 

® r 

sin ( ilf - e ) - 

- sin ( 71,-0 )j 

(3) 

r w 

r = V cos ( ijf - 6 ) “ cos 

(T1 -0)] 

(4) 


tan 


(5) 

o 

o 

It 

•> 



(6) 

h = 

c 



(7) 


where 

= magnitude of wind 

T1 = direction from which wind is blowing 
The three control inputs are: 

0^ = bank angle command 

a^ = longitudinal acceleration command 
h^ = vertical speed command 

3.3 SPIRAL ANGLE AS INDEPENDENT VARIABLE 

It will be convenient to replace time t with turn angle 0 as the independent 
variable in Equations (4) - (7). Dividing Equations (4) - (7) by (3) gives: 

^ = i= r cos - (W/V) cos (n-p) (S) 

d© % sin (^-©) - (W/V) sin (T] -©) 
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( 9 ) 


lt=l = L2. ° 

de e y2 sin (^-e) ^ (W/V) sin (l^-e) 


dv_ V _ ''^cA 

de S sin (^-9) - (WA') sin (n -9) 


( 10 ) 


rh /V 
c' 


h_^ 

d9 g sin (i;(-9) “ (WA) sin ('T]“9) 


( 11 ) 


The fifth state variable becomes time, which is defined by the reciprocal of Equation (3): 


dt _ 1 _ r/V 

dg 9 sin (i|f-9) - (W/V) sin (t] - 9) 

3.4 TIME IN SPIRAL FOR ARBITRARY TURN ANGLE 


( 12 ) 


One problem in specifying the nominal spiral and associated guidance 
parameters, is determining the nominal time as a function of turn angle, t*(0). The 
time to turn from 9^ to 9^ is (see Figure 9 for nomenclature): 



(1.3) 


where the ground speed V is found by the law of cosines: 

9 


(V^) = V 


Vl - cos^ (II - W sin (T1 - 0) 


(14) 


V V 

For constant spiral radius and airspeed, and defining the angle a = 9 - T| + tt/2: 
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(15) 



The elliptic Integral requires a series expansion to evaluate (Ref. 2): 


NORTH 



Figure 9 , Nomenclature for Determining t* (9), 
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E(-^/a) “ ® J ("v) f sin^0d0- |yA) J sin"^ ^ 


_ T3’5* -(2n-3) 
2*4*6* -211 

where the following relation can be used: 


(-^r i 


% . 2 n 

/ sin 


^ j _ sin^*^ ^orcosa l 2n-l f' 

0d 0 - 2^^ + y 


. 2 n r 
sin 0 cf 0 


. 2 n -2 ^ 
sin 0 d 0 


Thus, the elliptic integral can be written as 


q) “ sin O' cos O’) 

- (-^ 0 ! ~ -gl- sin a co« 


3 3 . 1.3 

(' O' “ ■ sin O' cos O' “ ' ^ 2 " " ® eos o?) 


+ H.O.T. 


For a complete turn, (i.e.),fi^ - 0 ^ = 2rr 




and the time for one complete turn is 


T = 2lli_ r,7w\2-', 1 /W\2 3 /J«\' 


W ° W 1 

Neglecting the ("^) f®rm in the series ex^i^sion of E (assuming — ^ < -j- 


causes an error less than 0 . 1 %; while neglecting the ( 77 ) term causes an error less than 

W 2 ^ 

0.3% . However, neglecting the (-rr-) term can cause an error greater than 6 %. 

V W 2 

Thus, it seems reasonable to keep terms out to (-^) and expect an accuracy of better 
than 0.3% . 

To verify this conclusion, the elliptic integral expansion was calculated for 
various values of (W/V) and o', using 2,3 and 25 terms. The results are compared in 
Table 3. Note that by keeping two terms in the expansion, i .e. (W/V) terms, 
the maximum error is only0.04%for W = V/2, i .e., sin ^(W/V) = 30°. 

- 20 - 


ABROBPACB ■YBTEIViai INC. • ONE vine brook park • burunston, masbacmusettb oiBoa * (ai7) a 7 e- 7 Bi 7 



Table. 3. Approximate Solution to Elliptic Integral of the Second Kind. 


Two Terms in Expansion 




0 

10 

20 

30 

40 

50 

60 

70 

80 

90 

0 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

10 

0.1745 

0.1745 

0.1744 

0,1743 

0.1741 

0.1740 

0.1738 

0.1737 

0.1736 

0.1736 

20 

0.3490 

0.3488 

0.3482 

0.3473 

0.3461 

0.3449 

0.3438 

0.3428 

0.3422 

0.3420 

30 

0.5235 

0.5229 

0.5209 

0.5178 

0.5140 

0,5100 

0.5061 

0.5029 

0.5008 

0.5000 

40 

0.6981 

0.6965 

0.6920 

0.6850 

0.6763 

0.6668 

0.6577 

0.6501 

0.6451 

0.6434 

50 

0.8726 

0.8697 

0.8614 

0.8483 

0.8318 

0.8138 

0.7964 

0.7818 

0.7721 

0.7688 

60 

1.0471 

1.0425 

1.0289 

1.0076 

0.9805 

0,9506 

0,9215 

0.8970 

0.8807 

0,8750 

70 

1.2217 

1.2149, 

1.1949 

1.1633 

1.1230 

1.0781 

1.0341 

0.9969 

0,9720 

0.9633 

80 

1.3962 

1.3869 

1.3597 

1.3164 

1.2607 

1.1985 

1.1371 

1.0850 

1.0500 

1.0376 

90 

1.5707 

1.5588 

1.5238 

1.4680 

1.3959 

1.3149 

1.2348 

1.1666 

1.1206 

1.1044 


Three Terms in Expansion 
0 sin ^(w/V) 



0 

10 

20 

30 

40 

50 

60 

70 

80 

90 

. 0 

O.OOQO 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

10 

0.1745 

0.1745 

0.1744 

0.1743 

0.1741 

0.1740 

0.1738 

0.1737 

0.1736 

0.1736 

20 

0.3490 

0.3488 

0.3482 

0.3473 

0.3461 

0.3449 

0.3438 

0.3428 

0.3422 

0*3420 

30 

0.5235 

0.5229 

0.5209 

0.5178 

0.5140 

0.5099 

0.5060 

0.5028 

0.5007 

0.5000 

40 

0.6981 

0.6965 

0.6920 

0.6850 

0.6762 

0.6667 

0.6575 

0.6498 

0.6447 

0.6429 

50 

0,8726 

0.8697 

0.8614 

0.8483 

0.8317 

0.8134 

0.7956 

0.7805 

0.7704 

0.7669 

60 

1.0471 

1.0425 

1.0289 

1.0075 

0.9802 

0.9496 

0.9193 

0.8934 

0.8759 

0.8697 

70 

1.2217 

1.2149 

1.1949 

1.1631 

1.1222 

1.0758 

1.0294 

0.9892 

0.9617 

0.9520 

80 

1.3962 

1.3869 

1.3596 

1.3161 

1.2593 

1.1944 

1.1286 

1.0712 

1.0317 

1.0176 

90 

1.5707 

1.5588 

1.5238 

1.4675 

1.3938 

1.3087 

1.2219 

1.1454 

1.0926 

1.0737 


Twenty-Five Terms in Expansion 


sin (w/V) 



0 

10 

20 

30 

40 

0 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

10 

0.1745 

0.1745 

0,1744 

0.1743 

0.1741 

20 

0,3490 

0.3488 

0.3482 

0.3473 

0.3461 

30 

0.5235 

0.5229 

0.5209 

0,5178 

0,5140 

40 

0.6981 

0.6965 

0.6920 

0.6850 

0.6762 

50 

0.8726 

0.8697 

0.8614 

0.8483 

0.8317 

60 

1.0471 

1.0425 

1.0289 

1.0075 

0.9801 

70 

1.2217 

1.2149 

1.1949 

1.1631 

1,1220 

80 

1.3962 

1.3869 

1.3596 

1.3160 

1.2589 

90 

1.5707 

1.5588 

1.5237 

1.4674 

1.3931 
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50 60 70 80 90 

0.0000 0.0000 0.0000 0,0000 0,0000 

0.1740 0,1738 0.1737 0.1736 0.1736 

0.3449 0.3438 0.3428 0,3422 0.3420 

0,5099 0.5060 0.5028 0.5007 0.5000 

0.6667 0.6574 0.6497 0,6445 0,6427 

0.8133 0.7953 0.7800 0,7697 0.7660 

0.9492 0.9183 0.8914 0.8727 0.8660 

1,0749 1.0266 0.9829 0,9514 0,9397 

1.1925 1,1224 1,0564 1.0057 0.9856 

1.3055 1.2110 1.1184 1.0421 1.0098 
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3.5 


NOMINAL SPIRAL SPECIFICATION 


The guidance concept in Figure 7 requires the specification of the nominal 
state vector x* and control vector u* as functions of the independent variable 6; i.e., 

x*^) = [r*(e), V*( 0 ), r(e)/ 1 *( 0 ), h*( 0 )]^ (22) 

u*{0) = [0*(0), o*(0), h*(0)]*‘' (23) 

The basic assumptions used for the nominal spiral are 

• Constant radius (r*) 

• Constant airspeed (V*) 

• Constant rate of descent (f\*) 

• Steady wind (T]*, W*) 

For right-hand turns (clockwise from above), the nominal spiral state 
variables are defined by the following equations: 


r* = 

constant 


(24) 

V* = 

constant 


(25) 

i*(e) 

= 0 +— - sin ^ cos (0 - T|*)] 

2 V* 


(26) 

t*(0) 


Q + '^ [cos (0 - 
^ V* 

T|*) 



+ cos(0 - T]*) sin (0 

-11*) 


- COS (0Q - T]*) sin (0 q - Tj*)] | 


(27) 

h*(0) 

= h*^Q) + h*[t*(e) "t*(0Q)j 


(28) 


where 0 q is the entry angle Into the spiral . 

The nominal heading. Equation (26), is obtained geometrically from 

Figure 8 . The nominal time expression uses the expansion of Equation (19), 

2 

neglecting terms of higher order than (W/V) and inserts this into Equation (15b). 


The nominal control variobles for a right-hand spiral are: 

1.2 r 


where 


0*(0) - tan 

A A [l 


-Ij V 


**• W* 

— A + IL sin (0 - Tl*) 

I r*gA L V* 


(29) 


\v*/ 


cos (0 - Tj*) 


1/2 
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a* = 0 


(30) 


h*(e.) - h*(0Q) 

h* = = constani- (31) 

t*(9f) - 


9^ is the exit angle from the spiral. The nominal bank angle is obtained by differen- 
tiating Equation (26) with respect to 0 , and substituting into Equation (5). 


Right Versus Left Turns 

For right turns, 9 “ 6 q ^ 0, and Equations (26) - (31) above apply. 
However, for left turns, 0 - 9 q < 0, and the following sign changes are required in 
Equations (26), (27) and (29): 


TT 


il»*(9) = 0 - — + sin 

2 


-1 ryr 

V* 


cos (0-11*) 


(26') 


3.6 


t*(9) = 


-il fi 

V*| \V*/ ( 


le-9o- 


w* 

V* 


[cos (0 - IT) -cos(0q-T1*)] 


- — [9 “ 9« + cos(0 - 1]*) sin (0 - 1]*) 

4\V*/ ^ 

- cos (0Q - T]*) sin (0 q - Tl*)]| 


0*(9) 


- tan 


1 v*^ 

A - ^ sin (0 - Tl*)l 1 

! r*gA 

V* J ' 


WIND EFFECTS ON NOMINAL SPIRAL 


(27-) 

(29') 


Except for r* and V*, all the nominal spiral parameters depend on the 
nominal wind speed W* and direction T]*. To illustrate the effects of wind on the 
nominal spiral. Table 4 compares several constant airspeed turns (V* = 60 kt, 
r* = 2000 ft) of 360°. The table shows the nominal time, ground speed, heading, 
heading rate and bank angle as functions of the spiral turn angle for winds from 0 
to 30 kt. Notice that in a 30 kt wind (W/V = 0.5), the time for the turn is more than 
30 seconds greater than the zero-wind turn. The ground speed varies from 30 kt to 
90 kt, the wind crab angle (i(t - 0) reaches 30°, and the heading rate exceeds 6.5 
deg/sec with a peak bank angle of nearly 20°. 
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Table 4. Wind Effects on 360° Turns (V —60 kt, R 2000 ft). 


Wind Speed = 0 kf 


9 

t 

Vg 



(deg) 

(sec) 

(kt) 

(deg) 

(deg/sec 

0.0 

0.0 

60.0 

0.0 

2.90 

30.0 

10.3 

60.0 

30.0 

2.90 

60.0 

20.6 

60.0 

60.0 

2.90 

90.0 

31.0 

60.0 

90.0 

2.90 

120.0 

41.3 

60.0 

120.0 

2.90 

150 0 

51.6 

60.0 

150.0 

2.90 

180.0 

62.0 

60.0 

180 .0 

2.90 

210 J3 

72.3 

60.0 

210.0 

2.90 

240.0 

82.7 

60.0 

240.0 

2.90 

270.0 

93 .0 

60.0 

270 0 

2.90 

300.0 

103.3 

60.0 

300.0 

2.90 

330 0 

IIB 7 

60.0 

330.0 

2.90 

360.0 

124.0 

60.0 

360.0 

2.90 


I 

ro 

4 ^ 

Wind Speed = 20 kt 
. . V. 


9 

(deg) 

r 

(sec) 

(kt) 

f w 

(deg) (deg/sec 

0.0 

0.0 

40.0 

0.0 

1.28 

30.0 

15.2 

41.8 

20.4 

1.43 

60.0 

29.2 

47.4 

43.2 

1.89 

90.0 

41.3 

56 .5 

70.5 

2.73 

120 0 

51.3 

67.4 

103.2 

B 82 

150.0 

59.9 

76.4 

140.4 

4.78 

180.0 

67.8 

79.9 

180.0 

5.15 

210.0 

75.6 

76.4 

219.5 

4.78 

240.0 

84.2 

67.4 

256.7 

3.82 

270.0 

94.3 

56.5 

289.4 

2.73 

300.0 

106.3 

47.4 

316.7 

1.89 

330.0 

120.3 

41.8 

339.5 

1.43 

360.0 

135.6 

40.0 

360.0 

1.28 


Wind Speed - 10 kt 


(deg) 

9 

(deg) 

t 

(sec) 

Vg 

(kt) 

(deg) 

(deg/ sec) 

0 

(deg) 

9.05 

0.0 

0.0 

50.0 

0.0 

2.01 

6.31 

9.05 

30.0 

12.3 

51.1 

25.2 

2.11 

6.62 

9.05 

60.0 

24.1 

54.3 

51.7 

2.40 

7.53 

9.05 

90.0 

35.0 

59.1 

80.4 

2.86 

8.93 

9.0 L 

120.0 

45.1 

64.3 

111.7 

3.37 

10.50 

9.05 

i5o:o 

54.4 

68.4 

145.2 

3.78 

11.76 

9.05 

180.0 

63.3 

69.9 

180.0 

3 94 

12.24 

9.05 

210.0 

72.2 

68.4 

214.7 

3 .78 

11.76 

9 .0 ! 

240.0 

81.6 

64.3 

248.2 

B 37 

10.50 

9.05 

270.0 

91.6 

59.1 

279.5 

2.86 

8.93 

9.05 

300.0 

102.6 

54.3 

308.2 

2.40 

7.53 

9.05 

330.0 

114.4 

51.1 

334.7 

2.11 

6.62 

9.05 

360.0 

126 7 

50.0 

360.0 

2.01 

6.31 


Wind Speed = 30 kt 


0 

(deg) 

9 

(deg) 

t 

(sec) 

(kt) 

(deg) 

(deg/sec) 

0 

(deg) 

4.05 

0 0 

0.0 

30.0 

0.0 

0.72 

2.28 

4.49 

30.0 

20.2 

32.1 

15.5 

0.85 

2.70 

5.94 

60.0 

37.9 

39.0 

34.3 

1.36 

4.29 

8.54 

90.0 

51.8 

51.9 

60 0 

2.51 

7.86 

11.88 

120.0 

62.1 

69.0 

94.3 

4.26 

13.19 

14.71 

150 0 

70.2 

84.0 

135.5 

5 88 

17.91 

15.82 

180.0 

77.3 

90.0 

180.0 

6.52 

19.73 

14.71 

210.0 

84.3 

84.0 

224.4 

5.88 

17.91 

11.88 

240.0 

92.4 

69.0 

265 6 

4.26 

13.19 

8.54 

270.0 

102.7 

51.9 

300.0 

2.51 

7.86 

5.94 

300 0 

116.6 

39.0 

325.6 

1.36 

4.29 

4.49 

330.0 

134.3 

32.1 

344.4 

0.85 

2.70 

4.05 

360 .0 

154.6 

30.0 

360.0 

0.72 

2.28 



If is apparenf from Table 4 thaf the nominal spiral is very sensitive to the 
predicted wind. If the wind speed and/or direction differ from the values used in the 
nominal calculations, the guidance system performance can suffer considerably. Con- 
sequently, a wind estimator will be an essential element in the guidance concept to 
update the forecast winds and to correct for wind shears during the descent. The 
estimated wind will then be used to compute the nominal bank angle and heading 
since they are the primary variables ilhat affect the desired radius. However, the 
nominal time will be based only on the original wind estimate to maintain proper air- 
craft spacing. The wind estimator is discussed in Section 5 of the report. 


3.7 FEEDBACK GUIDANCE LAWS 


As shown in Figure 7, the navigation system estimate of the vehicle 
state x(0) is compared with the nominal state x*(0), and any differences are fed back 
through the guidance gain matrix to generate corrections Au to the nominal control 
u*(0); i.e., 

u(0) = u*(e) + K [x{0) - x*(9)] (32) 


Theoretically, K is a 3 x 5 matrix so that an error in any element of the 
state will produce a correction in all three control channels. In practice, however, 
many of the elements of K can be neglected. As a result, the guidance laws of 
Equation (32) can be written for the individual control channels as follows: 


0^ = 0*(e) ± C^(r - r*) + C^U - r (0)1 + C^[t - t*(e)l (33) 

= K^(V-V*)+K^[t>-t*(e)l (34) 

he = h* + Kh [h-h*(0)l (35) 


The choice of upper or lower sign in Equation (33) depends on whether the turn is 
to the right or the left, respectively. 

For passenger comfort and safety, the guidance commands would be limited 
to reasonable levels. Typical constraints might be 


1 0, 

1 < 30° 

(36a) 

1 

1 < 0.1 g 

(36b) 


1 < 1000 ft/min 

(36c) 
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The selection of appropriate values for the feedback gains in Equations (33) - 
(35) is an iterative design process. A linearized analysis, described in Section 4, 
was used to provide a satisfactory set of guidance gains for preliminary evaluation 
of the guidance concept. 

The gains Cj, and Kj. are included in Equations (33) and (34) to provide 
4-D guidance capability. The guidance system can control time either with airspeed 
adjustments or by small changes in the radius, since the time varies with both; i.e. 



The change in final time is 



At V = 60 kt and r = 2000 ft, a I -knot change in airspeed has the same effect on final 
time as a 35 ft change in radius. Another consideration is that the airspeed response 
to commands is approximately linear with time, while the radius change is proportional 
to the second power of time (through ilf). Thus, time should be controlled only by 
airspeed adjustment and not by changing the radius. 


-26- 
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SECTION 4 

PERTURBATION GUIDANCE ANALYSIS 


4. 1 SIMPLIFIED 3-D GUIDANCE WITHOUT WIND 

A simplified analysis was conducted for 3-dimensional guidance during a 
zero-wind descent, with constant groundspeed (V = V ) and flight path angle. For 

w 

this simple spiral, the state variables are spiral radius r, spiral azimuth 0, aircraft 
heading and altitude h (Figure 10). The guidance scheme can be simplified by 
using 6 as the independent variable in place of time, and by separating the vertical 
and horizontal channels. 
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we can use 0 as fhe independent variable. Hence, dividing Equations (39) -■ (41) by 
(42), the perturbation equations of motion become: 


d(6ilf) 

= r*n^ 60 

(43) 

■ de 

c 


d(ftr) 

= -r* 6^ 

(44) 

d0 



d(6h) 

JA 

o 

• JZ 
•o 

II 

(45) 


where n^ ^ g (sec^ 0*)/V^ 

4.1.1 HORIZONTAL PLANE GUIDANCE 

The perturbation state and control vectors for motions in the horizontal plane 

are: 


X - 


6r 


, u = 60 




where x satisfies the linear differential equation 
^ = Fx + Gu 


where 


d0 ■ 


F = 


(47) 


0 0 


G = ’ 


r*n^ 


(48) 


Applying the Quadratic Synthesis technique, the optimum linear feedback 
guidance law is u = -Cx, or 


60 = -C . 6i)f - 6r 

* -IT 

where the feedback gain matrix isC = [C,C] = B G S 

■ i]; r 

The B matrix is determined by the maximum allowable bank control, ft0m: 


(49) 


(50) 


The steady-state gains are obtained from the solution of the matrix-Ricatti equations: 

S = -SF - F^S + SGB"^ G^S - A = 0 (51) 


where the matrix A describes the maximum allowable state perturbations: 
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A = 


(6t„) 


-2 


(«r ) 


-2 


m' 


The guidance gains then become: 


S = 




(rad/rad) 


(52) 


(53) 


6 0 

C = (rad/ft) 

r 6r 


(54) 


An alternate approach to selecting the feedback gains is provided by substituting 
the guidance equation into the perturbation equations of motion to obtain the characteristic 
equation of the closed-loop system: 


(6r)" +r*n^ C^, ( 6r)' + r*Vc ( «r ) = 0; 

where 

(Sr). = -L 

r* d9 


(55) 


(56) 


This is a second-order system whose angular natural frequency and damping ratio are 
given by 




(57) 


* 2^ 
2tc 


(58) 


n 

A reasonable response can be obtained by specifying the frequency of response in terms 
of a characteristic angle, 0^^ of the order of 1/4 of the turn, and a damping ratio ^of 


approximately 0.7: 
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(59) 


= (2tt/ 9^) /n^r*2 




(60) 


Thus, any radial error would be correcfed wifh critical damping after a turn angle 

of approximately 9 . 

c 

4.1.2 VERTICAL PLANE GUIDANCE 

This is a scalar problem, with the state and control defined from Equation 
(45) as follows: 


X = 6h; u = 6 h 


( 61 ) 


F = 0; G = r*/V 


(62) 


The feedback guidance law is 


= -'<h 


6h 


(63) 


where the feedback gain is 


K = (bh )^ S 
h ' c'm Y 


(64) 


The Ricatti equation is a scalar: 


S = (r*A)^ 0 

cm m 


(65) 


Thus the feedback gain is 

(6h ) 

\ r*' I 


Kh = 


c' m 


6h 


(1/sec) 


m 


( 66 ) 
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Using the alternate approach, the characteristic equation is 


(6h)' +I1 k. (6h) = 0 
V ^ 

Choose the characteristic response angle here to be 9 ; i.e., 

Kh = 

4.1.3 EXAMPLE 

As an example of the feedback gains required, assume a 
following parameters: 

r* = 2000 ft 
V* = 60kt = I01.27fps 
V = -6° 
i = 3°/sec 

Selecting the characteristic angle to be 6^ = tt/2, give 

K = 9 . 443 “ 
r? = 0.003224 ft"' 

Lfct 


60 = 
m 

5° 

( 6h ) 

c-'m 

= 200 fpm 

5r = 

75 ft 

6h 

= 100 ft 

m 

m 


- 

m 

15° 




These give the guidance gains: 

Cj. = 0.06667 deg/ft 
= 1.0272 rad/rad 

K|^ = 0.03333 ft/sec/ft 
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Using the alternate approach/ the same example gives: 


= 0.07111 deg/ft 
C ^ = 0.8685 rad/rad 
= 0.03224 ft/sec/ft 

Comparison of these results indicotes very good agreement for the two approaches. 

4.2 4-D GUIDANCE 

The perturbation state equations for no wind, linearized about the nominal 
spiral are: 

d ( 6 x)/d9 = F &x + G 8 u (69) 



The feedback guidance laws of Equations (33) ~ (35) are, in vector-matrix 

form: 

5u = -C 6x (71) 

where 

o Ky o K}. o 

-C=Cr o C\(i C|. o 

_o o o o 

A block diagram of the horizontal guidance system is shown in Figure 1 1 ., 
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Figure 11 . Linearized Block Diagram of Horizonfal Guidance System. 


The characteristic equation for this system is given by 

det (Ip - F + G C) = 0 (73) 

In expanded form this gives 



For constant airspeed, Ky = Kj = V = 0, and the characteristic equation for 
the horizontal axes becomes 



AEROaPACE aYBTEIVIB, INC. • ONE VINE BROOK PARK • BURUNSTON, MASBACHUBETTB aiBOS • (B17) a7S-7B17 



For s^ability, fhe guidance gains must satisfy the following constraints 


C / C. > 0, and C , < 0, and C. < ^ C 

Approximate magnitudes of and were determined in Section 4. 1 .3 for r* = 2000 ft 
and V* = 60 kt: 

= 1 rad/rad 

C =0.001 rad/ft 
r 

This requires C^< 0.635 rad/sec. To achieve approximately 0.7 damping, we can select 
= 0.1 rad/sec 

Now it is convenient to convert back to the time domain: s A d/dt = B d/d© = 
(V*/f*) P* The above gains give 

s^ + 0.3177s +0.03217S+ 0.001609 = 0 = (s^ + 2CU)j^ s + u) (s+ 1 /t) (76) 

which yields tw = 0.092 rad/sec, ^ = 0.66, t — 5. 6 sec. 

The performance of the bank guidance loop with no acceleration feedback and 
perfect navigation was examined by simulation for the first two minutes of a level turn*. 
Table 5 shows the nominal time histories, while Tables 6-8 illustrate the effects of initial 
errors in r, 'jf and V. Note that nominal conditions have been re-established within one 
half a turn for the first two cases; however, since there is no acceleration feedback, the 
airspeed error creates steady-state radius and time errors. Tables 9 and 10 present the 
effects of a 5-knot wind from 90® and 180°; it is worth noting that even such small winds 
cause significant radial position errors. Table 11 shows that a 20-knot wind from the 
north produces very large errors in position and time. However, these errors can be 
almost eliminated, as shown in Table 12, by including the wind in the nominal calculations. 


*Section 8 describes the simulation programs. 
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Table 5. Nominal Turn With Bank Control Only 


t 

e 

'If 

R 

0.0 

0.0 

0.0 

2000.0 

5.0 

14.5 

14.5 

2000.0 

10.0 

29.0 

29.0 

2000.0 

15.0 

43.5 

43.5 

2000.0 

20.0 

58.0 

58.0 

2000.0 

25.0 

72.5 

72.5 

2000.0 

30,0 

87.0 

87.0 

2000.0 

35.0 

101.5 

101.5 

2000.0 

40.0 

116.0 

116.0 

2000.0 

45.0 

130.5 

130.5 

2000.0 

50.0 

145.0 

145.0 

2000.0 

55,0 

159.5 

159.5 

2000.0 

60.0 

174.0 

174.0 

2000.0 

65.0 

188.5 

188.5 

2000.0 

70.0 

203.0 

203.0 

2000.0 

75.0 

217.5 

217.5 

2000.0 

80.0 

232.0 

232.0 

2000.0 

85.0 

246.5 

246.5 

2000.0 

90.0 

261.1 

261.1 

2000.0 

95.0 

275.6 

275.6 

2000.0 

100.0 

290.1 

290.1 

2000,0 

105.0 

304.6 

304.6 

2000.0 

110.0 

319.1 

319.1 

2000.0 

115.0 

333.6 

333.6 

2000.0 

120.0 

348.1 

348.1 

2000.0 



f * 



9.05 

0.00 

0.0 

9.05 

9.05 

5.00 

14,5 

9.05 

9.05 

10.00 

29.0 

9.05 

9.05 

15.00 

43.5 

9.05 

9.05 

20.00 

58.0 

9.05 

9.05 

25.00 

72.5 

9.05 

9.05 

30,00 

87.0 

9.05 

9.05 

35.00 

101.5 

9.05 

9.05 

40.00 

116.0 

9.05 

9.05 

45.00 

130.5 

9.05 

9.05 

50.00 

145.0 

9.05 

9.05 

55.00 

159.5 

9.05 

9.05 

60.00 

174.0 

9.05 

9.05 

65.00 

188.5 

9.05 

9.05 

70.00 

203,0 

9.05 

9.05 

75.00 

217.5 

9.05 

9.05 

80.00 

232.0 

9.05 

9.05 

85.00 

246.5 

9.05 

9.05 

90.00 

261.1 

9.05 

9.05 

95.00 

275.6 

9.05 

9.05 

100.00 

290.1 

9.05 

9.05 

105.00 

304.6 

9.05 

9.05 

110.00 

319.1 

9.05 

9.05 

115.00 

333.6 

9.05 

9.05 

120.00 

348.1 

9,05 


R* = 2000 ft 
V* = 60 kt 
W* = Okt 
n* = Odeg 

= 0,001 rad/sec 
C, = -1.0 
Cj. = 0 rad/sec 
R^ = 2000 ft 
= Odeg 
V = 60 kt 

o 

= 0 deg 
W = 0kt 
T| = 0 deg 


Table 6 . 


Effect of Initial Radial Error with Bank Control Only. 


t 

e 

Ilf 

R 

0.0 

0.0 

0,0 

2200,0 

10.0 

26.7 

36.3 

2052.8 

20.0 

55.7 

58.4 

1940.8 

30,0 

85.7 

84.5 

1932.1 

40.0 

115.6 

113.9 

1960.6 

50.0 

144.9 

143.9 

1985.2 

60.0 

174.1 

173.7 

1997.6 

70.0 

203.1 

203.0 

2001.4 

80.0 

232.1 

232.1 

2001.6 

90.0 

261.1 

261.1 

2000.9 

100.0 

290.1 

290.1 

2000.3 

110.0 

319.1 

319.1 

2000.0 

120.0 

348.1 

348.1 

1999.9 


d > 

t * 


0 * 

20.51 

0.00 

0.0 

9.05 

6.89 

9.22 

26.7 

9.05 

7.46 

19.21 

55.7 

9.05 

8.83 

29.56 

85.7 

9.05 

9.36 

39.84 

115.6 

9.05 

9.34 

49,97 

144.9 

9.05 

9.20 

60.01 

174,1 

9.05 

9.10 

70,01 

203.1 

9.05 

9.05 

80,01 

232.1 

9o05 

9.04 

90.00 

261.1 

9.05 

9.04 

100.00 

290,1 

9.05 

9.05 

109.99 

319.1 

9.05 

9.05 

119.99 

348.1 

9.05 


Rq = 2200 ft 
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Table 7, Effects of Initial Heading Error With Bank Control Only 


t 

' 6 


R 


t* 

f 

0* 

0.0 

0.0 

10.0 

2000.0 

-0.94 

0.00 

0.0 

9.05 

10.0 

29.4 

27.1 

1971,2 

8.90 

10.13 

29.4 

9.05 

20.0 

58.5 

57.3 

2005.9 

9.54 

20.18 

58.5 

9.05 

30.0 

87.3 

87.4 

2015.1 

9.19 

30.12 

87.3 

9.05 

AO.O 

116,1 

116.5 

2010.9 

9.01 

40.05 

116.1 

9,05 

50.0 

145.0 

145.3 

2004.9 

8.98 

50.01 

145.0 

9.05 

60.0 

174.0 

174.1 

2001.3 

9.01 

59.9 

174.0 

9,05 

70.0 

203.0 

203.1 

1999.9 

9.03 

69.99 

203,0 

9.05 

80.0 

232.0 

232.0 

1999.6 

9.05 

79.99 

232.0 

9.05 

90.0 

261.0 

261.0 

1999.7 

9.05 

89. 9^ 

261.0 

9.05 

100.0 

290.1 

290.1 

1999.8 

9.05 

99.99 

290.1 

9.05 

110.0 

319.1 

319.1 

1999.9 

9.05 

110.00 

319.1 

9.05 

120.0 

348.1 

348.1 

2000.0 

9.05 

120.00 

348.1 

9.05 



Table 8. Effects of Airspeed Error With Bank Control 

Only 

t 

9 

♦ 

R 

<t> 

t* 



0,0 

0.0 

0.0 

2000.0 

9.05 

0.00 

0.0 

9.05 

10.0 

31.2 

27.9 

2036.1 

10.05 

10.76 

31.2 

9.05 

20.0 

61.5 

58.7 

2099.2 

10,55 

21.22 

61.5 

9.05 

30.0 

91.2 

89.7 

2140.6 

10.31 

31.43 

91.2 

9.05 

40.0 

120.4 

119.9 

2158.7 

10,02 

41.50 

120.4 

9.05 

50.0 

149.4 

149.3 

2164.6 

9.87 

51.52 

149.4 

9.05 

60.0 

178.5 

178.4 

2166,1 

9.81 

61.53 

178.5 

9.05 

70.0 

207.5 

207.5 

2166.5 

9.79 

71.53 

207.5 

9.05 

80.0 

236.5 

236.5 

2166.5 

9.79 

81.53 

236,5 

9.05 

90.0 

265.5 

265.5 

2166.6 

9.79 

91.53 

265.5 

9.05 

100.0 

294.5 

294.5 

2166.6 

9.79 

101.53 

294.5 

9.05 

110.0 

323.5 

323.5 

2166.6 

9.79 

111,53 

323.5 

9.05 

120.0 

352.5 

352.5 

2166.6 

9.79 

121.53 

352.5 

9.05 
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Table 9 . Effect of 5-Knot Wind With Bonk Control Only (90°). 


t 

0 


R 


t* 

i|i* 

6 * 


0.0 

0.0 

0.0 

2000,0 

9.05 

0.00 

0.0 

9.05 


10.0 

27.9 

32.1 

2053.7 

10.01 

9.64 

27,9 

9.05 


20.0 

54.7 

62.2 

2010.1 

8.67 

18.86 

54.7 

9.05 


30.0 

81.7 

88.3 

1911.6 

7.77 

28.18 

81.7 

9.05 


AO.O 

110.2 

112.7 

1817.6 

7.57 

37.98 

110.2 

9.05 

W = 5 kt 

50.0 

140.4 

137.5 

1770.7 

8.00 

48.39 

140.4 

9.05 

60.0 

171.8 

164,3 

1787.7 

8.81 

59.23 

171.8 

9.05 

T] = 90 deg 

70.0 

203.6 

193.8 

1860.7 

9.66 

70.21 

203.6 

9.05 

80.0 

235.1 

225.9 

1967.3 

10.29 

81.07 

235.1 

9.05 


90.0 

265.8 

259.5 

2079.8 

10.56 

91.63 

265.8 

9.05 


100.0 

295.2 

293.4 

2169.6 

10.44 

101.78 

295.2 

9.05 


110.0 

323.4 

326.3 

2214.5 

9.98 

111.47 

323.4 

9.05 


120.0 

350.3 

357.3 

2203.8 

9.32 

120.77 

350.3 

9.05 




Tab! 

e 10. 

Effect of 5- 

-Knot Wind With 

Bank Control 

Only (180°). 

t 

0 


R 

0 

t* 

\tr* 

6 * 


0.0 

0.0 

0.0 

2000.0 

9.05 

0.00 

0.0 

9.05 


10.0 

31.0 

28.9 

2044.5 

9.64 

10.71 

31.0 

9.05 


20.0 

60.5 

61.3 

2123.3 

10.31 

20.87 

60.5 

9.05 


30.0 

88.2 

93.7 

2149.0 

9.71 

30.41 

88.2 

9.05 


40.0 

114.7 

123.4 

2103.7 

8.79 

39.56 

114.7 

9.05 


50.0 

141.1 

150.3 

2008.7 

8.02 

48.65 

141.1 

9.05 

W = 5 kt 

60.0 

168.3 

175.2 

1898.0 

7.60 

58.01 

168.3 

9.05 

T) = 1 80 deg 

70.0 

196.9 

199.4 

1807.2 

7.60 

67.88 

196.9 

9.05 

80.0 

227.1 

224.3 

1764.4 

8.05 

78.31 

227.1 

9o05 


90.0 

258.6 

251.2 

1781.3 

8.81 

89.14 

258.6 

9.05 


100.0 

290.4 

280.7 

1851.6 

9.61 

100.11 

290.4 

9.05 


110.0 

322.0 

312.6 

1956.3 

10.24 

110.99 

322.0 

9.05 


120.0 

352.7 

346.0 

2069.2 

10.55 

121.59 

352.7 

9.05 
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Table 

11. Effechof20- 

Knot Wind With Bank Control Only. 



f 

e 


R 


t* 


(t>* 



0.0 

0.0 

0.0 

2000.0 

9.05 

0.00 

0.0 

9.05 



10.0 

19.7 

30.5 

1844.8 

7.63 

6.82 

19.7 

9.05 



20.0 

44.0 

48.5 

1508.3 

4.00 

15.18 

44.0 

9.05 



30.0 

78.4 

60.6 

1315.4 

4.55 

27.04 

78.4 

9.05 



40.0 

117.9 

82.3 

1455.8 

9.64 

40.66 

117.9 

9.05 



50.0 

155.2 

120.5 

1828.1 

13.78 

53.50 

155.2 

9.05 

W* 

= 0kt 

60.0 

189.3 

168.3 

2268.9 

15.26 

65.28 

189.3 

9.05 

W 

= 20 kt 

70.0 

219.6 

217.1 

2620.8 

14.40 

75.69 

219.6 

9.05 

80.0 

245.1 

260.5 

2774.5 

12.19 

84.50 

245.1 

9.05 



90.0 

266.2 

296.1 

2710.9 

9.78 

91.76 

266.2 

9.05 



100.0 

283.7 

324.1 

2462.9 

7.11 

97.80 

283.7 

9.05 



110.0 

298.9 

346.1 

2074.1 

5.99 

103.03 

298.9 

9.05 



120.0 

313.5 

362.8 

1586.7 

4.39 

108.08 

313.5 

9.05 





Table 12. 

Effect 

of Including Wind Estimate in 

Nominal 

Calculation. 

t 

0.0 

0 

0.0 

0.0 

R 

2000.0 

4.05 

t* 

0,00 

0.0 

(ff* 

4.05 


10.0 

19.4 

13.0 

19 9.9 

4.24 

9.99 

13.0 

4,23 


20.0 

39.6 

27.4 

1999.5 

4.86 

19.98 

27.3 

4.83 


30.0 

61.5 

44.5 

1998.6 

6.07 

29.97 

44.4 

6.04 


40.0 

86,2 

66.9 

1997.6 

8.20 

39.96 

66.8 

8.16 


50.0 

115.3 

97.8 

19 6.9 

11.42 

49.94 

97.8 

11.36 

W* = 20 kt 

60.0 

70.0 

149.8 

188.0 

140.2 

190.7 

1995.9 

1995.1 

14.76 

15.72 

59.93 

69.95 

140.1 

190.7 

14.70 

15.73 

W = 20 kt 

80.0 

225.3 

238.9 

1996.2 

13.34 

80.00 

239.0 

13.41 


90.0 

257.7 

276.4 

2000.5 

9.80 

90.07 

276.8 

9.86 


100.0 

285.0 

303.3 

2007.0 

7.09 

100.14 

303.7 

7.10 


110.0 

308.3 

323.3 

2012.6 

5.46 

110.15 

323.5 

5.42 


120.0 

329.3 

339.2 

2013.7 

4.55 

120.11 

339.1 

4.51 
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Turning now fo the airspeed control loop, it is apparent from Figure 11 that 
it decouples from the bank control loop if Cj, = 0. In this case, the characteristic 
equation of the acceleration loop, in the time domain, is: 

s^ - Kys + K/V* = 0 (77) 

Thus the natural frequency is with damping C ~ "Kv/2^n* Selecting the 

gains to provide approximately the same response as the bank loop gives 

Ky = -0.15 sec ^ 

Kj. = l.Oft/sec^ 

Using these acceleration loop gains and the bank loop gainspreviously defined, the 
response to an airspeed error is very acceptable, as shown in Table 13. The response to the 
same error, but with = 0 is shown in Table 14; note that the radial position error is 
smeller than before, while the time error is about the same. The responses to a 20-knot 
wind for C^=0. 1 and = 0 are compared in Tables 15 and 16 . On the average, the 
timing errors are about the same, but the radial error is considerably improved by the 
decoupled system. However, the radial error still exceeds 300 ft, which would generally 
not be acceptable. 

It is readily apparent from this analysis that the winds are a very important 
disturbarice on the spiral guidance system performance. In order to reduce the sensitiv- 
ity to wind effects, it will be necessary to include the wind estimates in the nominal 
path calculations, and perhaps increase the guidance gains. 
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Table 13* Effects of Airspeed Error With Coupled Bank & Throttle Control. 


f 

e 


R 

V 


a 

t* 



0.0 

0.0 

0.0 

2000.0 

50.00 

9.05 

1.500 

0.00 

0.00 

9.05 

10.0 

27.0 

29.6 

1959.4 

59.65 

8.09 

0.461 

9.30 

27.00 

9.05 

20.0 

57.2 

55.9 

1947.2 

61.38 

8.87 

-0.048 

19.73 

57.24 

9.05 

30.0 

87.2 

85.2 

1982.2 

60.38 

9.57 

-0.103 

30.07 

87.25 

9.05 

40.0 

116.3 

115.7 

2005.6 

59.74 

9.36 

-0.019 

40.09 

116.33 

9.05 

50.0 

145.1 

145.3 

2006.7 

59.79 

9.05 

0.021 

50.01 

145.10 

9.05 

60.0 

174.0 

174.2 

2001.0 

59.99 

8.97 

0.013 

59.97 

174.00 

9.05 


K^= 0.15 sec ' 

K = 1 .0 ft/sec^ 
' -1 
C^ = 0.1 sec 

V =50kt 
o 




Table 14. Effects of Airspeed Error With Uncoupled Bank & Throttle Control. 


t 

0 


R 

V 


a 

t* 


0 * 


0.0 

0.0 

0.0 

2000.0 

50.00 

9.05 

1.500 

0.00 

0.00 

9.05 


10.0 

26.9 

26.8 

1985.6 

59.68 

8.32 

0.474 

9.27 

26.92 

9.05 


20.0 

56.5 

55.7 

1998.1 

61.82 

9.77 

0.017 

19.50 

56.59 

9.05 


30.0 

86.3 

86.2 

2007.2 

61.39 

9.61 

-0.073 

29.77 

86.37 

9.05 

C^ = 0.0 sec 

40.0 

115.7 

115.9 

2006.4 

60.71 

9.26 

-0.054 

39.91 

115.78 

9.05 

V = 50 kt 

50.0 

144.9 

145.1 

2003.4 

60.30 

9.10 

-0.027 

49.96 

144.96 

9.05 

o 

60.0 

174.0 

174.1 

2001.5 

60.11 

9.06 

-0.011 

59.99 

174.03 

9.05 




Table 15. Effects of Wind With Coupled Bank & Throttle Control. 


t 

0 


R 

V 

d> 

a 

t* 




0.0 

0.0 

0.0 

2000.0 

60.00 

9.05 

0.000 

0.00 

0.00 

9.05 


10.0 

20.8 

29.6 

1846.4 

65,72 

7.62 

0.817 

7.17 

20.80 

9.05 


20.0 

50.3 

45.5 

1576.4 

70.70 . 

4.66 

-0.047 

17.37 

50.39 

9.05 


30.0 

87.1 

62.8 

1568.9 

65.49 

8.47 

-0.847 

30.03 

87.14 

9.05 


40.0 

122.0 

93.2 

1725.9 

57.00 

10.23 

-0,786 

42.08 

122.09 

9.05 


50.0 

155.1 

129.8 

1923.9 

50.10 

10.09 

-0,567 

53.46 

155.10 

9,05 

C^ = 0.1 sec 

60.0 

185.5 

170.6 

2154.3 

45.90 

10.13 

-0.235 

63.96 

185.57 

9.05 

o 

<N 

II 

70.0 

212.7 

213.5 

2365.6 

45.51 

10.05 

0.187 

73.35 

212.79 

9.05 

80.0 

236.8 

254.2 

2479.1 

49.42 

9.70 

0.606 

81.65 

236.88 

9.05 


90.0 

258.7 

288.3 

2431.2 

57.10 

8.81 

0.913 

89.19 

258.75 

9.05 


100.0 

280.4 

313.6 

2214.7 

66.75 

7.20 

0.958 

96.67 

280.45 

9,05 


110.0 

305.5 

331.0 

1920.8 

74.82 

5.80 

0.553 

105.31 

305.52 

9.05 


120.0 

336.0 

346.4 1729.8 77.13 

Table 16^^. Effects of Wind 

7.07 -0.103 115.83 336.05 9.05 

With Uncoupled Bank & Throttle Control. 


t 

0 


R 

V 


a 

t* 




0.0 

0.0 

0.0 

2000.0 

60.00 

9.05 

0.000 

0.00 

0,00 

9.05 


10.0 

20.7 

19.9 

1924.7 

65.72 

5.53 

0.829 

7.14 

20.74 

9.05 


20.0 

47.6 

38.7 

1837.7 

72.08 

8.63 

0.313 

16.41 

47.61 

9.05 


30.0 

80.6 

66.0 

1793.3 

71.62 

11.79 

-0.435 

27.79 

80.62 

9,05 


40.0 

117,8 

101.4 

1781.7 

64.27 

12.96 

-1.018 

40.63 

117.89 

9.05 


50.0 

156.1 

142.9 

1818.4 

52.84 

11.79 

-1.182 

53.80 

156.10 

9.05 

C^ = 0,0 sec 

60.0 

190.5 

185.6 

1913.6 

42.56 

9.08 

-0.759 

65.69 

190.58 

9.05 

W= 20 kf 

70.0 

218,9 

224.0 

2055.2 

38.48 

7.10 

-0.004 

75.45 

218.90 

9.05 

80.0 

242.3 

256.5 

2202.3 

41,71 

6.42 

0.641 

83.54 

242.37 

9.05 


90.0 

263.5 

283.3 

2297.0 

50.05 

6.24 

0.992 

90,84 

263.54 

9.05 


100.0 

284.5 

305.0 

2309.0 

60.49 

6.24 

1.055 

98.09 

284.57 

9.05 


110.0 

307.1 

323.9 

2255.4 

70.42 

6.79 

0.891 

105.85 

307.10 

9.05 


120.0 

331.9 

342.8 

2171.9 

78.00 

7,99 

0.598 

114.43 

331.97 

9.05 
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4.3 


GUIDANCE GAINS 


The linearized model of the horizontal guidance law in Figure 11 allows the 
analyst to apply all the standard techniques of linear control theory to select an 
appropriate set of guidance gains. However, the following cautions should be noted: 

• The actual guidance low places a limit on the magnitude of the 
commanded bank ongle. 

• The preceding analysis has neglected any attitude dynamics in the 
rotorcraft's response to bank angle commands. 

• The actual system may operate outside the range where the linear 
model is valid. 

To avoid potential problems, the linear model is used to select the guidance gains and 
to predict system performance; and the prediction is then verified using the non-linear 
model in a trial and error fashion. Improved gains were selected using this technique 
together with a determination of the effect of a finite bank angle rate. 

The analysis led to selection of the following guidance gains: 

C^. = 0.005 rad/ft 

= -1.77 rad/rad 
= 0 rad/sec 

The value of is based on results presented in Section 3.7; has a relatively 
minor effect on the 4-D guidance performance, and tends to destabilize the 3-D response. 
The guidance gains C^, and were selected iteratively from combinations that provide 
a natural damping of 70%. Figure 1 2 compares the radius error, heading error, and 
bank histories in response to a 100 ft initial radius error for four sets of gains. The error 
responses of the first three cases are very similar, while case four is extremely slow in 
settling out. Case 3 has a much less violent bank history, and was therefore selected as 
the best compromise. 

In Figure 1 2 the actual bank angle begins at the maximum allowable value 
(30°) for cases 1 - 3 in response to the initial condition error on r. To determine the 
effect of a finite bank angle rate, the guidance low was modified to generate a bank 
rate command proportional to the error between desired and actual angle, viz 
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The results for = 10 sec ^ are compared in Figure 1 3 with the infinite bank rate 

results of Figure 12 . The finite bank rate has little effect on the radial error history/ 
and produces only a small delay in the heading response. For lower values of the 
response lags more, but is practically the same within 40° of turn angle for = 1,0 sec ^ 


4.4 


ERROR SENSITIVITY 


Figure 1 1 is valuable for determining the sensitivity and dynamic response of 
the system to measurement errors. Disturbances can be entered at the points labelled 
6 0, 6Y, 6r, 6V, and ^t to correspond, respectively, to errors in the measurement of 
bank angle, heading angle, radial position, airspeed and time. Since the steady-state 
input to all the integrators has to be zero, the position error sensitivity to heading 
error, for example, is given by 


£r 


steady state 



= 6 ft/deg 


(79) 


= .3 deg/ft 

= 1.8 deg/deg 


1 ^ 

cu 

n 


2.5 sec 


^ = 7.5 deg 5 = .7 

P 


(80) 


Similarly, the steady-state sensitivity of position error to bank angle error is given by 

= 3.5 ft/deg. (81) 

60 Cr 
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SECTION 5 
WIND ESTIMATION 


The previous section has demonstrated the need for a reliable estimate of the 
winds in order to achieve the required guidance accuracy during the spiral descent. 
This section presents two approaches to providing this capability. 

5. 1 KALMAN WIND ESTIMATOR FORMULATION 

The first approach is an extension of the wind estimator in the Kalman filter 
developed for the Straw-Man navigation system (Reference 1). The nomenclature for 
the Kalman wind estimator for spiral approaches Is shown in Fig e 14. 



Figure 14. Nomenclature for Wind Analysis. 

The North and East wind components are 
Wj^ = -W cos 

We ~ “W sin Tl 

The tangential and radial components are: 

# 
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(83a) 


W^. = -W sin (T) - 0) = sin 0 + cos 9 

= -W cos (T) - 0) = W|^ cos 0 + Wg sin e (83b) 

The state equations (8) - (12), in terms of the North and East wind components. 


are 


dr _ 


— r 


do 


V cos (t - 0) + W|^ cos 0 + sin 9 

V sin (\|f - 0) - sin 0 + cos 0 

tan 0 


^ = rg 

do V V sin (tj/ - 0 ) - VVj^ sin 0 + cos 0 


dV _ 


ra. 


de 

V sin (\|f - 

0) - W|,^ sin 0 + cos 0 

dh ^ 


rh 

c 

d0 

V sin (ijf - 

0) - sin 9 + Wg cos e 

dt ^ 


r 


d0 


V sin (t ” 6) “ W|>^ sin 0 + cos 0 


(84) 

(85) 

( 86 ) 

(87) 

( 88 ) 


These must be augmented by the equations for the wind component models. For example, 
if these are assumed to be exponentially-correlated random processes ( i.e., first-order 
colored noises), the appropriate equations are 


d9 


(89a) 

dWn 

We 


U 

d0 

- — + q£ 

0_ 

(89b) 


where 0^ is the characteristic angle of the processes, and q^ and q^ are driving white 
noises. 

The optimum estimator for the wind Is formulated os follows. The linearized 
state equations with the guidance loop closed are: 
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dx _ c , 
— = Fx + q 

de 


(90) 


The perturbation state vector is 


X = [6r, 6ijr, 6V, 6h, 6t, 


(91) 


where 6W|^ = North component of wind error 
6Wg = East component of wind ^rror 

The driving noise vector is 


q = [0, 0, 0, 0, 0, q^, q^I^ 


(92) 


The linearized closed-loop system matrix is approximately 


0 

-r 

0 

0 

0 

r cos s/V 

+r sin g/V 

rgC/V^ 


0 

rgC/V^ 

0 

+rg tan 0 sin 0/V^ 

2 

-rg tan 0 cos 0/V 

0 

0 

rK/V 

rK/V 

0 

0 

0 

0 

0 

0 

0 

rK^A 

0 

0 

1/v 

0 

-r/V^ 

0 

0 

+r sin e/V^ 

-r cos g/V^ 

0 

0 

0 

0 

0 

-1/Qc 

0 

0 

0 

0 

0 

0 

0 

-1/e^ 


The linearized measurement vector is given by 


m = Hx + r (94) 


where r is the vector of measurement noises, and H is the measurement geometry matrix. 
For example, if direct measurements are made of 6r, 6V, 6^, Stand 6h, then 



0 

1 

0 


0 0 

I 

0 i 

I 

1 ! 

1 ! 

I 

1 ! 

0 0 

0 0 


(95) 
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The estimate of the state x is propagated by 


X = Fx + PH^R“^ (m - Hx) (96) 

where R is the power spectral density of r. The error covariance matrix P propagates as 
fol lows; 

P = FP + PF^ + Q - PH^R“’ HP (97) 

where Q, the power spectral density of q, defines the wind strengths. 

Further development and analysis of the Kalman wind estimator were discontinued 
in favor of the simplified wind estimator described subsequently. 


5.2 


SIMPLIFIED WIND ESTIMATOR 


One simplified wind estimator can be formulated by reasoning that the lateral 
errors 6r and 6\|f are due primarily to the radial wind component, while the longitudinal 
errors 6V and fit are due primarily to the tangential wind component: 


(98) 


z 

to 

1 


sin 6 cos 0 


' K^fir + 

K2 

I 

LU 

1 


-cos 0 sin 0 


+ 

1 

et_ 


The net effect of this estimator is to provide an additional integration to null the "quasi 
steady-state" error due to an unknown wind. However, the proper gains are difficult 
to pick by trial and error, so it may still be easier to use the optimum Kalman formulation. 

Another simplified wind estimator is obtained by using the off-nominal values 
of bank angle and heading angle to estimate the tangential and radial wind components. 
The changes in the north and east components are transformed from the tangential and 
radial components by the spiral angle 0. 


1 

o» 


-sin 0 

cos 0 


o> 

I 

1 

LLJ 

<o 

1 


cos 0 

sin 0_ 


1 

u. 

to 

1 


If the rotorcraft Is close to the nominal spiral, a tangential wind requires a 
change in bank angle, while a radial wind requires a change in heading into the wind. 
Thus, the bank and heading perturbations give an indication of the unknown wind: 


- 50 - 


AEROSPACE SYSTEMS, INC. • one VINE brook park • Burlington, mabbachusettb oisoo • (bit) s7e-7Bi7 


( 100 ) 


69 \ B0 / 

0 V / 


6W 

60 


6 0 


6\)f 


(101) 


where 0^ is a "characteristic turn angle" over which the wind estimate is averaged, and 
the estimator "gains" are given by 

aw. 


30 


^ 

2V* 


"bank sensitivity to tail wind" 


( 102 ) 


aw 

L = V* "crab sensitivity to crosswind" (103) 

aiji 

Example 

r* = 2000 ft 


V* = 60 kt 

For a characteristic time of the estimator of about 30 sec, let 0 = 90®. The estimator 

c 

gains are 

I!9 ~ 320 - 3kt/deg (104a) 

2V* rad 

V*— 1 00 — ~ 1 kt/deg (104b) 

rad 


The performance of the simplified estimator is shown in Figures 15, 16 and 17 for 
an unknown steady wind of 20 kts from four directions, using the above gains and the 
same nominal spiral as used previously. The wind speed estimate is within 62% of the 
actual after a 90® turn in all cases, and the direction estimate Is even better. The 
maximum radial error is less than 90 ft for initial cross winds, and less than 50 ft for an 
initial head or tail wind. 

From these results, the simplified estimator seems to provide an effective means 
of reducing the effects of unknown winds on the spiral guidance scheme. 


- 51 - 


AEnOSPACE SYETEMEi INC. • one vine BROOK PARK • BURLINSTON, MASBACHUBETTB OIBOS • (B17) S7S-7S17 


ESTIMATED WIND SPEED, W (KT) 



Figure 15. Estimated Wind Speed With Simplified Estimator, 
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ESTIMATED WIND DIRECTION, 7 ) (DEG) 



Figure 16. Esfimafed Wind Direction With Simplified Estimator. 
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240 270 300 330 


_J 

360 


ied Wind Estimator 


5.3 WIND ESTIMATOR DYNAMICS 

Because knowledge of the wind is crucial to the performance of the guidance 
law, it is important to consider the dynamics of the wind estimator. The effect of the 
wind estimator is to produce variations in 0* and i)r* to eliminate radial errors: 

_ 1 


PV 


(♦ - D / ^0 " ° 


(105a) 


p0 


* = 


e 


{<P^~ 0*) / ~ constant 


(105b) 


where p = d/d6. Substituting (105a) and (105b) into (33) with = , solving for 0 

gives 


0 = C / 1 + 
'x r 




+ C.I, t 


(106) 


where 6r = r - r*. 


The result of the wind estimator is to produce an effective integral -plus-bypass 
term in the feedback of 6r. This is exactly the technique used in mechanizing ILS 
couplers for conventional aircraft. 

The linear signal flow diagram in the "angle domain" is shown in Figure 18. 
The characteristic equation of this system is third order; 


+ Cp + c/*^9 ’ 




Y*'" V*^ ®< 


For stability, the characteristic angle of the estimator must satisfy 

.2 


0 c >- 


V’ 


r*gC, 


(107) 


(108) 


The sensitivity of the estimator to a crosswind is 

"C, w 

6r S 

C V* 
r 

For Cj. = 0.005, = -1.77, V* = 60 kt, this gives 6r = 6 ft/kt of unknown 

crosswind. 


(109) 
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Figure 18. Linear Signal Flow Diagram for Wind Estimator. 

The wind estimator will lag in a steady shear wind by approximately the 
characteristic angle 0^. The maximum error due to a shear is about 

w 0 hX_(dW/dh) fir (110) 

c 

Conventional automatic landing systems are normally certificated on the basis of dW/dh = 
8 knots/30 meters « 0.08 kt/ft '^'With 0^ = 30° and li = 500 ft/min, this produces a 
maximum radial error in a shear of 38 ft. Figure 19 illustrates the performance of the 
wind estimator during a 300° spiral descent of 1000 ft in a steady shear of 0.01 kt/ft. 

As expected, the estimated wind lags the true wind, with a peak error of about 3 kts. 

For straight flight during a transition segment (see Section 7), the preceding 
analysis is modified as follows: 

0^= CyX + C^ (♦-)*) (111) 

A 

t* = —(?-'!'*)= (112) 

9cP 1 

+This is not maintained over altitude changes of more than a couple hundred of feet. 
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WIND SPEED, KT 


Figure 19. Wind Estimator in a Steady Shear 

Therefore, 



QqP 

1 +e^P 




The linearized block diagram is shown in Figure 20. 


(113) 


- 57 - 


AEROSPACB aVBTEMB, INC. • ONE VINE BROOK PARK • BURLINBTON, MASSACHUSETTS C1B03 • (817) S7S-7B17 


Figure 20. Linear Signal Flow Diagram for Wind Estimafor During 
Transition Segment. 
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SECTION 6 

PHASE PLANE ANALYSIS 


In this section the nonlinear response of the guidance laws are Interpreted 
in the phase plane. The problem of capture when spiral guidance Is initiated a long way 
from equilibirum conditions leads too modification in the guidance algorithm. 

6.1 THE PHASE PLANE 

A * 

The phase plane is a plot of heading error, t , vs radial position error, 
r - r*. Since the lateral guidance variable, 0^, is determined as a nonlinear function of 
these two parameters, it is possible to plot lines of constant bank in the phase plane as 
shown in Figure 21 . It is also possible to plot trajectories of the nonlinear response of the 
system as shown in Figure 22. 
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Some examples of horizontal system performance for large initial errors are shown in 
Figure 23. In examples A and B, large initial heading and radial errors, respectively, 
are eventually corrected and the rotor craft trajectory settles into the desired spiral in 
an equilibrium condition. In example C, the extremely large initial radial error pre- 
vents the rotorcraft from leaving the maximum right-turn bank command, and consequently 
it fails to capture the equilibrium spiral condition. 


6.2 GUARANTEED SPIRAL CAPTURE 

For Initial spiral offsets greater than 1200 - 1500 ft, the spiral guidance laws 
will command a continuous maximum bank angle (left or right) and result in a circular 
flight which never intercepts the desired spiral. Although such large deviations should 
not normally occur, the guidance law must still guarantee capture in the event they do 
exist. 

6.2.1 RIGHT HAND TURNS 


To ensure spiral capture, the guidance law can be modified to restrict the 
heading angle error, . The modified phase plane is shown In Figure 24 for 

driving Uf - \lr MS A \li . The modified bank guidance algorithm becomes 

a i T I a a 




,r>r, 




0* + C^(r - r*) + C^(? - 1 ^*) ,r 2 <r<r^ 


C. (i - t ) 

V max 


r < r 

, ' Ir 


where 


C I A i +0 

r = r* T 

ri " 1 


r2 = r* + 


C , A t - 0' 
^ max 


and is the largest value of 1 I’ - 'll*! for 0^ = 0. 


(114) 


(115) 

(116) 
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Figure 23 . Phase-Plane Examples. 
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Figure 25 shows two examples of the modified guidance law performance with 

a large initial error for =60° and 90°. The capture region in Figure 25 is only 

shown for A =±60°. There is a different capture zone for the case of A| =+90°. 
max max 

The 60° case has a large initial heading error overshoot of hut relatively smaller 

final overshoots of the nominal 1 fading and radius. The initial heading overshoot for 
the case of Ai = +90° is smaller. If the capture region were shown it would be 
evident that the initial overshoot does not exceed the linear region. The final overshoots 
are larg«er in the case of ~ +90°. Consequently, 60° seems to provide a good 

compromise solution . 

6.2.2 LEFT-HAND TURNS 

For left-handed spirals, the phase plane is effectively "flipped" around the 

A 

- '1'* axis. The modified bank guidance law is 
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^c = 


c , “ 'll + A ijr ) 


,r>r, 


0* - C^(r - r*) + - 'll*) 


C I ('ll - 'll* - A 'll ) 
' 'll ^ max' 


/ r <r. 


where 



C,A'^ - 0' 

-il ^ max 



+ 


C I A'li + 0' 
'll max 


(114') 

(IIS') 

(116') 
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SECTION 7 

SPIRAL TRANSITION SEGMENTS 


The transitions between the spiral and straight flight are critical phases which 
must be accomplished smoothly. In general, the nominal transition paths will consist 
of straight line segments tangent to the spiral arcs. On the circular arcs, the inde- 
pendent variable is the spiral angle ( 9 ) measured around the center of the arc. Along 
the transition segments, the independent variable is the along-track distance (x) 
measured in the direction of the nominal ground track (Figure 26). 



Figure 26. Independent Variables for Straight Segments and Circular Arcs. 

7. 1 NOMINAL PATH FOR ENTRY TRANSITION TO SPIRAL 

7.1.1 RIGHT-HAND SPIRAL ENTRY 

A constant groundspeed, constant-altitude transition is assumed, which is 
tangent to the descent spiral . Slightly more general equations are developed for the 
simulation described In Section 8 for the straight segments. There, groundspeed Is 
assumed to vary linearly with time during periods of constant acceleration rather than 
remaining constant as is assumed in this section. An x-y coordinate frame is defined 
with its origin at the spiral center and with the x-axis parallel to, and in the same 
direction as, the transition path (Figure 27). The guidance logic switches at 0 = 9 q: 

9 < 6q — guide to straight transition segment 

9 ^ 9q — engage spiral guidance logic (^17) 
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Figure 27. Entry Transition to Right-Hand Spiral. 

For a constant groundspeed transition, the nominal path is defined as follows: 

V* = const. = [V*^ + W*^ + 2V* W* sin ( 11 * - (118) 

g g ' O'-" 

‘ff = const. = 9q + — sin ^ [(W*/V*) cos (Tl* - 0 q)] (1 19) 


y* = const. - -r* (120) 

t* = t*(x) = t*(6o) - VV * (Note, x<0) (121) 

0* = o (122) 

h* = 0 (123) 

h = const. = h* (0q) (124) 
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( 125 ) 


Equation (118) can be solved for in terms of V : 

V* = [V*^ + W* ^ - 2V*W* cos (Tl * - D] 
9 




7.1.2 LEFT-HAND SPIRAL ENTRY 

The x-y axes are defined in the same way as for a right-hand spiral, and the 
logic is also switched at S = 0 q (Figure 28); 

6 > — guide to straight transition segment 

(117') 

0 < 0Q — engage spiral guidance logic 

N 



Figure 28. Entry Transition to Left-Hand Spiral. 


For a constant-speed transition. Equations (118), (119) and (120) are modified as follows: 

V* = const. = [V*^ + W*^ - 2V*W* sin (Tl* - 0q)] (118') 

9 9 ^ 

r = const. = 0Q -j+sin"M(WVV*) cos (H* - 0q)] (119') 

y* = const. = r* (120') 
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7.2 NOMINAL PATH FOR EXIT TRANSITION FROM SPIRAL 
7.2.1 RIGHT-HAND SPIRAL EXIT 


This is essentially the same as the transition onto the spiral, with a constant 
deceleration to hover. As shown in Figure 29, a set of x-y coordinates is defined, 
centered at the nominal spiral axis with x parallel to the final approach transition 
segment. The guidance logic switches at 6 = 0|r: 

0 < 9^ — guide to spiral path 

0 > 0|T — guide to straight transition segment (126) 

N 



Figure 29. Exit Transition From Spiral. 

For this transition, the nominal airspeed and heading will not be constant. For a 
constant deceleration, the nominal path equations are: 
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where 


where 


t*(x) = tf- - t*(9j.)] V 1 - x/xj. 
+f‘'+*(e^) 


x^ = V*( 9 f) 


2 


(127) 


V*(x) = V*(9f) [ 1 - 


t*(x) - t*(9^) 


tf- t*(ep 


] 


( 15 a) 


V*(9f) = [V*^ + W*^ - 2V*W* cos (7]* - I 


9= 0r 


V*(x) = [V*(x)^ + W*^ + 2V*(x)W* sin (T]* - 9^)] 

(129) 

1 V*(x)A^* + sln (T]*-9r) 

ili*(x) = 9r + tan [ ^ ] 

cos (T]* - 9j:) 

(130) 

-] ( W*V*(9r) cos {T]* - 9r) ) 

0*(x) = - tan < § > 

( gV*(x)[t* - t*(9|:)] ) 

(13il) 

y* = -r* = const. 

(132) 

-V*(9f) 

a* = 9 [V*(x) + W* sin (T]* - 9r)] 

V*(x)[ff*-t*(9f)] 9 * 

(133) 

h* - h*(er) 

h* = = const. 

tf - t*(0f) 


h*(x) = hf + h* [t*(x) - tf] 

(135) 


- 71 - 


AEROSPACE BYBTEIVISi INC. • one VINE BROOK PARK • burunstdn, Massachusetts oiBos • (ai7) S7E-7Bn7 



7.2.2 LEFt-HAND SPIRAL EXIT 

The only differences between the left and right-hand spiral exits are the 
sign of 0, and 180 in the ground track angle at 0p The guidance logic switches at 9^: 

9> 0|; — guide to nominal spiral 



0 < 0 ^ — guide to straight transition segment 

(136) 

The modifications in the above equations are: 



V*(x) = [V*(x)^ + W*^ - 2V*(x)W* sin ('ll* - 0 r)] 
9 9* 

<T37) 


sin ( 11 * - 6f) - V*(x)A\I 

,|f*(x) = 9c - tan I — ] 

cos (T|* - ep 

(138) 


(W*V*( 0 n) cos (Tl* - 0 n)) 

0 *(x) = tan < ^ ) 

( gV*(x)[t* - t*( 0 p] ) 

(139) 


y* = r* = const. 

(140) 


-V*(9c) 

a* = i— 1 [V*(x) - W* sin ('ll* - 0 r)] 

V*(x)[tf- t*( 9 ,)I 9 ' 

(141) 

7.3 

TRANSITION SEGMENT GUIDANCE LAW 


and (35), 

For. the circular arc, the proposed guidance law is given by Equations (33), (34) 
For consistency, a similar guidance law is proposed for the straight segments: 


= + C^Cy - y*l + - ♦* (i)l 

(142) 


[V - V* W] + ft - t*Wl 

(143) 


= h*(x) + K|^[h - h*(J)l 

(144) 
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Note that, unlike Equation (3-32), the longitudinal acceleration command in Equation 
(143) contains a nominal term in addition to the feedback terms. 


Outside the linear bank command region (0=*f'0 ), the guidance logic must 

mox ^ ^ 

be modified to ensure capture of the transition segment. For large offset errors (y - y ), 
the logic is adjusted to provide a 60° intercept of the nominal path (? - t*); as shown in 
the phase plane of Figure 30, The change in the guidance algorithm from that of 
Equation (142) is: 


0 


C 


C^,[t-+*(x)-n/3] ;<y, 

0 (x) + Cy [y - y*l + C| ll - y^ S Jsy^ 

It- ♦* (x) +T1/3I y>y 2 


(145) 


where y ^ ~ y* 



- * a. 

X2-y + 




C 

y 
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For the entry transition the linear region is symmetric about the origin since 
0* = 0. This is nearly true in the exit transition as well since the nominal bank angle 
will be very small . The slope of the trajectories in the phase plane is 


^ _ y _ V*^ sin (t - 
^ g tan 0^ 


(146) 


Outside the linear region, 0^ - —^max' paths are given by 

- v*2 


(p _ y*) = + V cos (t - \ ) + (y - y*) 1 J ^ 


g tan 0 


90 ' 


(147) 


max 


7.4 HOVER AND TOUCHDOWN PHASE 

The final phase of the approach and landing is not unique to the spiral descent. 
A hover mode must be engaged when the rotorcraft is over the pad. Using a pad-centered 
frame (Figure 31), 


N 

4 


PAD 




♦ E 


Figure 31 . Pad-centered Frame for Hover and Touchdown. 


the hover mode is engaged when the rotorcraft is within a desired radius r^^ of the hover 
point; i .e . , when 


[N^ + E^] < r, « r 

n p 


(148) 


where r is the radius of the pad. The nominal conditions are 
P 
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f =T1* (149) 

V* = W* (150) 

^* = ^hover 051) 

h* = 0 (152) 

0* = O (153) 

a* = 0 (154) 


After a short hover (~5 sec), a slow vertical descent to the pad is performed 



= const. 


(155) 
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SECTION 8 
SIMULATION 


Two digital simulation programs have been developed for the analysis and 
evaluation of the spiral guidance concept. The firsts containing a simple model of the 
vehicle dynamics and no navigation errors, was used In an interactive mode to verify 
the algorithms and to investigate guidance gain adjustments. The second is a modifica- 
tion of the VALT simulation program described in Reference 1 . 


8.1 INTERACTIVE SIMULATION 


A relatively simple simulation, written in BASIC, was developed for ASI's 
inhouse minicomputer to provide preliminary verification of the guidance concept. 
This interactive simulation was particularly valuable for allowing the analysts to 
conveniently adjust the guidance algorithms, parameters and gains, and thereby gain 
additional insight and understanding of the system. 


This simulation numerically integrates the nonlinear equations of motion for 
the point-mass rotorcraft model. Equations (3) through (7); the finite bank angle rate 
command of Equation (78); and the following simplified wind estimator equations, from 
Equations (100) and (101): 



(156) 


057 ) 


The following limits are observed in the simulation: 


|0C 

I < 0 « 

'^max 

CO 

o 

o 

(158) 

|0| 

• 

<0 ^ 
^max 

lO^/sec 

(159) 

Ih 

c 

1 < h 

' max 

1000 ft/min 

(160) 

The nominal 

path varlabl 

es are calculated as functions of 9 

using Equations 


(24) through (31). The guidance commands are generated using the laws of Equations 
(34) and (35), with the modified bank command of Equation (117). 
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Since the simplified wind estimator is based on the premise of small perturba- 
tions in and 0 from their nominal values, the wind estimates are not updated if the 
rotorcraft is outside an elliptical region in the center of the phase-plane: 



(161) 


where 90° and Ar 200 ft. The estimator is also deactivated whenever the system 
is outside the linear bank command region of the phase plane. 


The interactive simulation can also be used to analyze the spiral transition 
segments discussed in Section 7. Thus, it was used to verify the algorithms for left and 
right spirals, for entry and exit transitions, for nominal time, for wind estimation, for 
guaranteed capture of the nominal, etc., as well as to select a set of satisfactory guidance 
gains. The resulting guidance algorithms were then implemented in the Program VALT 
simulation. 


8.2 PROGRAM VALT SIMULATION 

The second simulation is a modification of Program VALT (VTOL Automatic 
l-anding Technology) which was developed under a previous investigation (Ref. 1) to 
analyze the rotorcraft navigation system performance and to conduct subsystem sensitivity 
studies. The program is written in Fortran IV for operation on the LaRC CDC 6400/6600 
computer system. Figure 32 presents a general flov/ diagram of VALT to illustrate the 
overall organization and logical operation of the simulation. 

In essence, VALT is actually a double simulation. First, it integrates the 
equations of motion which describe the response of the rotorcraft and flight control system 
to the guidance system commands, and it simulates the actual noisy output of the various 
navigation sensors. This part of VALT is a nonlinear, stochastic process which is intended 
to provide a reasonably accurate representation of the "real world." The second part of 
VALT simulates the operations of the onboard navigation and guidance systems. This 
portion also models the rotorcraft motions and navigation measurements, but here the 
models are much simpler, and are linearized about a desired nominal flight path. The 
models used in this part of the program are purposely kept as simple as possible to minimize 
the onboard computation requirements. 

The principal outputs of the simulation are time histories of the rotorcraft 
position and velocity, and two sets of error histories. The estimator errors are the dif- 
ferences between the estimated and actual position/velocity, and thus indicate the 
navigation systems' performance. The second set of errors show the rotorcraft' s actual 
deviations from the nominal position/velocity profiles; these illustrate the overall 
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Ftgure-32. Overall Flow Diagram of Program VALT 
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performance of the entire system, including the navigation, guidance scheme, flight 
control system, rotorcraft capability, and wind effects. 

Program VALT was modified for the present investigation by appending to 
the existing guidance algorithms the option of a spiral entry transition, spiral descent and 
exit transition. This required the addition of the nominal path calculations, the wind 
estimator and the new guidance commands. The spiral guidance calculations are divided 
into three segments, as shown in Figure 33. 


NORTH 



Figure 33. Segments of Nominal Spiral. 

• Segment A-B: Entry Transition 

• Segment B-C: Spiral Descent 

• Segment C-D: Exit Transition 

A rectangular reference frame in the horizontal plane, centered at the pad is used (x,y). 
Figure 34 presents an overall flow diagram of the spiral guidance algorithm. Details of 
the guidance calculotions are summarized in the following subsections. 

8.2.1 INITIALIZATION CALCULATIONS 

The initialization phase is executed once at the beginning of the flight to 
calculate a number of parameters at the points indicated in Figure 33: 

• Point A is the last waypoint of the enroute path. The given navigation 
coordinates are its latitude, longitude, altitude, time, groundspeed and 
ground track: 
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INITIALIZATION 



NO 


YES 


FLAG= 1 


ISEG=I 

J 


FLAG = 0 

CALCULATE 9 q , 6 ^ , ETC 
I FOR RIGHT TURN 
I FOR LEFT TURN 




ENROUTE GUID. 
CALCULATIONS 



ENTRY TRANSITION 
ALGORITHM 


SPIRAL DESCENT 
ALGORITHM 


EXIT TRANSITION 
ALGORITHM 


HOVER 
CALCULATIONS 


CHD 

Figure 34. Overall Flow Diagram of Spiral Guidance Algorithm. 
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( 162 ) 


*-A' 

The guidance coordinates are: 

’‘A = - ‘‘dK ‘-D 

>'a “ ^‘•A ■ '■D^e 

where r is the earth's radius, 
e 

• Point B is the spiral entry point: The initial spiral angle is*: 


^A "^E 

, -1 

+ cos 

r* 

yA-ye 

,V(>‘E-''A>-^''<yE->'A>^ . 




(163a) 

(163b) 


(164) 


where 9 q > 0 for right turns and 6 q ^ for left turns. The guidance coordinates, 
assuming a constant deceleration along segment A-B, are: 


Xg = -i- r* sin 6 q 

yg = yg + r* cos 9g 


where 


l-B = 'a + 2V(x5-x^)^ + (y5-y//(Vg^ + 

V = [v*^' - W*- cos^ (e^ - 11*)] + w* sin (e^ - n*) 

gg u - u 


• Point C is the spiral exit point. The exit spiral angle is: 


6 ^ - tan ^ 


-Xr 


■yeJ 


X -1 

+ cos 


L>/^J 


(165) 

(166) 

(167) 

(168) 

(169) 


(170) 


where 0|; > 9 q for right turns and 9|^< 9 q for left turns. The guidance coordinates are: 
x^ = Xg + r* sin 6^ (171) 


Upper algebraic sign is used for right turns; lower sign for left turns. 
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+ r* cos 9^ 

he = hg + hg (te - fg) 




_2(w 

4\V* 


V* 

* \ 2 


1 - 


f W* 


-1 




(172) 

(173) 

T1*)-cos^q-T1*)] 

(174) 


[Qf-Qo"^ cos (0f-'n*)sin(9f-Tl*)-cos(9Q-Tl*)sin(eQ-Tl*)] 


The nominal rate of descent is assumed constant from point B to point D, I.e., 

If |hgj > spiral is increased by one turn {9 ^ = 0^ ±2tt), and Eqns. (170) to 

(176) are recomputed. 

m Point D is the hover point over the pad. The given navigation coordinates are: 
^D' ^D' ^D' ^gp ® 

The guidance coordinates^ assuming a constant deceleration from C, are: 

x^ = yD = 0 (176) 




J 2 , ' 2 

2\x^ + 


V 


9C 


V = [V*^ - W*^ cos^ (9r-n*)]^'^^ + W* sin (9p-T]*) 

a>- f — T 


(177) 


(178) 


Sc 


Point E is the spiral center relative to the pad, with guidance coordinates: 
^E' 


8.2.2 CALCULATING NOMINAL PARAMETERS 

The spiral guidance algorithm must compute the nominal path variables based 
on the estimated position x, y, fi-. 

• Entry Transition Segment (A) - (B). The estimated position relative to the 
spiral center is: j 


Ax = X - Xg 

(179) 

Ay = y - y^ 

(180) 


- 83 = 


AEROSPACE SYBTEMSi INC. • one vine brook park • Burlington, Massachusetts 01B03 • (bit) a 7 a-?Bi 7 


It is convenient to use a rotated coordinate system with x parallel to the entry 
transition segment and origin at the spiral center (Figure 35): 


X = £ (Ax cos 0Q - Ay sin 9 q) 

(181) 

y = T (Ax sin 0Q + Ay cos 9 q) 

(182) 


N 



N 



RIGHT TURN 


Figure 35. Rotated Coordinate Frame for Spiral Entry. 

If X > 0, the guidance logic switches to the spiral mode [Segment (B) - (C)] . The 
nominal guidance parameters are calculated along the entry transition as functions of x: 


t*(x) = 




tn + 


r / 

I -V 
^B 


\p- -V^ )-V ,V 

n_ 'n n.' r\ ^ r\ • i 


^B ^A 


— , = V 

V % 9A 

9A 




9b 9a 

(183) 


where = -y/(xg - Xa)^ + (yg - Yjf 


where 


y* = + r* 


V*(x) = lV*(x)^ + W^+ 2V*(3 i)Wsin(fi«-Ti )]^/2 

9 9 ^ 


V*(x) 

g 


(V -V ) 

9b 9a 

fL[tg-t*(x)] 


9 b (to - t . ) 


B 'A' 
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(188) 


il;*(x) == tan 


0*(x) = - tan I + 


_1 I ± cos 0Q + [ W/V*(x)l sin Tt 
+ sin 9 q + [WA*(x)1 COST] 

I W COS (6q “ T] ) 


gV*(x) 




N -V 

a*(x) = +1 — ;-) sin [9 q - t*(x)l 


'b“'a / 


h* = 0 




( 189 ) 


( 190 ) 

( 191 ) 

( 192 ) 


• Spiral Descent Segment (B) - (C) 

/ 

During the spiral phase, the nominal guidance parameters are determined as 

^ A 

functions of the estimated angle around the spiral, 9 . If 0>9p the guidance logic 
switches to the exit transition phase. The nominal parameters during the spiral are: 

( 193 ) 

( 194 ) 


r* = const. 
V* = const. 


TC 

t*(e) = tn± — 
D V* 




[cos (9 “ i]*) 


- cos (9 q -T]*)] - [e -9 q + cos (9 -T]*)sin(9 -n*) 


-cos(0q-T]*) sin(0Q-Tl*)] 


'' -1 
0*(9) - + tan 


A 

W * A - 

A ± sin (9 - T| ) 

2 


r*gA 

V* 



( 195 ) 


( 196 ) 
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where 


A = 1 1 


■(^) 


cos (0 ” 7] ) 


lV2 


a* = 0 


A A 


\ TT . -1 

j ~ - Sin 

yi COS (0 -T])l 

(2 

* 

> 


* 

h* = T r — = consf. 

X "’’B 


• A 


h*(e) = hg + h*[t*(e) - tg] 


( 197 ) 

098) 

(199) 

(200) 


• Spiral Exit Segment (C) - (D) 

A rotated coordinate system, centered at the pad is used for the exit segment 
(Figure 36): 

X = ± (x cos 6j; “ y sin 0j.) (201) 

y = +(x sin0|; + y cos 0^;) (202) 



Figure 36. Rotated Coordinate Frame for Spiral Exit. 
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V ~2 ~2 

X + y < r , the guidance logic switches to the hover modeo Otherwise the 
P 

nominal parameters on the exit transition are; 


t*(x) = t^ + (t^ - tj^) 


( 203 ) 


where 


~ - J 2 , 2 
’'c - * >'c 


y* = 0 


V*(x) = lV*(x)^+W^T2V (x)W sin (9f- 11)1 


(204) 

(205) 

(206) 


where 


V*(x) = V 

g 9 C 


tf. -t*(x) 


’’D ”^C 


(207) 


A A 


ijf*(x) = tan 


-1 


I ± V* (x ) cos 9|- + W sin T) 


+V* ( X ) sin 9|: + W cos f| 


(208) 


0*(x) = -tan 


-1 


W V cos (9n - T] ) 
gV*(x)(tj;,-tc) 


(209) 


a*(x) = ± 


V sin [9p - i|f*(x)l 
9r ' 




(210) 


h* = hp + h* [t*(x) “ tp] 


(211) 


■ ~ 

h* = = const. 


( 212 ) 
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• Hover Mode 

Hover guidance is not unique to the spiral approach, and can be analyzed 
separately. For the present study, the simulation terminates at the completion of the 
spiral exit transition, 

8.2.3 WIND ESTIMATES 


The simplified wind estimator described in Section 5 is engaged during Segment 
(B) - (C), when the aircraft is close to the nominal path: i.e.. 




< 1 


(213) 


which defines an ellipse about the nominal path in the phase plane. The wind estimator 

is also deactivated whenever I 0,- I > 0 

'^max 

When the wind estimator is active, changes in the radial and tangential wind 

A 

components are estimated over an increment in G : 

6W 


r _ V* 


60 

6W 




c 


&e 2v*e 


(214) 


(215) 


These are converted into navigation coordinates: 



A 

6W^ 

6W, 


COS 0 

- si n 9 

t 

(216) 

89 


60 

80 


A 


A 

A 


6Wp 

A 

6W 

6W, 


E _ 

sin 9 

+ COS 9 

t 

(217) 

89 


60 

60 
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The North and East wind estimates are then updated: 


A A 


w ^(9 + 60 ) = (e) + 


A A 


Wg(e + 69 ) = W^( 9 )+ I ^ j 69 


(•-?) 

n) 


69 


The estimated magnitude and direction of the wind are obtained from: 

w = Vw^ + 

A 

T] = tan 

8.2.4 GUIDANCE LAWS 



(218) 


(219) 


( 220 ) 


(221) 


The spiral guidance commands are summarized below for the three flight 

segments: 

• Spiral Transition Segments (A) - (B) and (C) - (D) 


t C^[1'-^*(x)“A'lr^Q^] ,y<(y*-Ay) 

0*(x) + Cy[y - y*l + C^[ J-\!(*(x)], (y* -Ay)< y<(y* + Ay) (222) 

C^U-r(x) + A ,y>(y*+Ay) 



a^ = a*(x)+Ky[V-V*(x)]+K^[t-t*(x)] 

(223) 


he = li* + h - h* (x)] 

(224) 

where 

Ay - Ai|f C,/C 

' ’^max Y y 

(225) 
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Spiral Descent Segment (B) - (C) 

,r>r, 

A A A A A 

= |0»{e)±C^(r-r*) + C^[4-t*(9)l, r,a r2rj 


c^Ui-r(e) + A*„^^l 


/ r< f 2 


( 226 ) 


where 


°c " K^[V-V*l + K^[t-t*(9)l 
h = h* + K. [h- h*(e)l 

c n 




C, All; ± 

ax 


0*(e) 


(227) 

(228) 
(229) 


C.Ailr +0*(e) 
= r* + ^ 


(230) 


• Guidance Command Limits 


I 0 I < 0 « 30° 

’^c ’^max 


I a I < a 2 kt/sec 

c max 

I h I ^ li w 1000 ft/min 
c max ' 


(231) 

(232) 

(233) 


• Guidance Gains 

Suggested values for the feedback guidance gains are summarized in Table 17. 
8.2.5 MEASUREMENTS 

The guidance system inputs required for the previous calculatiof ,5 are assumed 
to come from the navigation system and other onboard equipment. 
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Table 17. Recommended Values 
for Guidance Gains. 



0.005 

rad/ft 

s 

-1.77 

rad/rad 

Cf 

0 

rad/sec 

K 

-0.15 

-1 

sec 

V 

S 

1.0 

ft/sec^ 

Kh 

0.033 ft/sec/ft 


• Position 

A ^ A 

Position estimates are obtained from the navigation system outputs of L, Jl, h: 



(234) 

A A 

y = (L-Lp)r^ 

(235) 

" „ Jr s2^ " v2 

r - ^(x - Xg) + (y “ yg) 

(236) 

A - A A 

e = tan [(x - Xg)/(y-y^)] 

(237) 


• Airspeed and Heading 


Airspeed and heading estimates are obtained from the navigation systems 

A A A 

inertial velocity estimates (Vj^, V^, Vp) and the wind estimator: 



- W^) + (Vf - W,:)^ + 


Pa a 



(238) 

(239) 


- 91 - 


AEROSPACE SYSTEMSi INC. • one VINE BROOK pahk • Burlington, Massachusetts oibq 3 • (aiv) a 7 a- 7 si 7 




• Attitude 


Attitude estimates are obtained in the simulation from the lateral acceleration: 


0 = tan 


-1 


cos t sin f 


mg 


(240) 


where Fj^ and F^ are the North and East components of total force , 


• Time 

Time estimates are provided by a clock with a bias A tg : 

t = t + Atg (241) 

where the time bias is selected from a zero-mean, Gaussian distribution with a 
standard deviation a. « 2 sec. 

8.3 SPIRAL APPROACH EXAMPLE 

To check out the spiral guidance algorithms in Program VALT, and to evaluate 
their performance, a realistic example has been chosen for the New York City Kennedy 
International Airport (JFK)o Figure 37 shows two example spiral approaches to an 
existing helipod (point D) at JFK, The entry transition (A - B) overflies the helipod at 
1 100 ft; a 270° descending spiral is made to the right (B - C) or left (B - C*); and a 
short exit transition (C - D or C - D) completes the maneuver to hover (D). 
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Figure 37. Example Spiral Approaches to Kennedy International Airport (JFK) 




The VALT simulation with the spiral guidance system has been used to de- 
termine the effect of a variety of initial condition errors, wind errors, guidance 
parameter variations and navigation system errors for the example approach to JFK. 

The nominal initial conditions for the runs are: 

L = 40.7053 N 

o 

I = 73.7803 W 

o 

h = 1100 ft 

o 

V = 120 kt 

go 

= 180 deg 
= 186.6 deg 
W = 15 kt 

T| = 90 deg 

These correspond to a point 2.67 nm north of point A, the beginning of the entry 
transition segm.ent. 

Without navigation measurement errors, the spiral guidance system follows the 
nominal path almost perfectly. The maximum horizontal error is about 150 ft and occurs 
four seconds after the entry transition segment logic is engaged. It is corrected within 
12 seconds to less than five feet, and remains negligible throughout the descent. This 
is merely a turn-on transient effect which is the result of a discontinuity in the nominal 
ground track when switching from the enroute guidance mode to the spiral guida/ice 
algorithm. The position and velocity error histories are depicted in Figures 38 and 39. 

The remaining VALT runs investigated various errors and system parameter 
variations. These are summarized in Tables 18 & 19. Typical results are depicted in 
Figures 40 through 48. In general, large initial condition errors can introduce signi- 
ficant errors due to nonlinear effects. 
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0 < 
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Table 18, Spiral Guidance Parametric Analyses. 


INITIAL CONDITION ERRORS 

Cross -track 

+ 1 nm 

Along-track 

+ 1 nm 

Ground speed 

+ 20 kt 

Ground track 

+ 30 deg 

Altitude 

+ 300 ft 

WIND ERRORS 

Wind speed (constant) 

+ 15 kt 

Wind direction (constant 

+ 90 deg 

Wind speed (shear) 

+ 1.36 kt/100 ft 

Wind direction (shear) 

+ 8.18 deg /l 00 ft 

GUIDANCE PARAMETERS 

Radial error gain {Cp 

0.001, 0.010 rad/ft 

Heading error gain (C^) 
Airspeed error gain (K^) 

-1.0, -3,0 rad/rad 
-0.10, -0.30 sec“^ 

Altitude error gain (K^) 

0,5, 2 ,0 sec“^ 

Characteristic angle (0^) 

45°, 180° 

Nominal wind speed (W*) 

+ 15 kt 

Nominal wind direction (T]*) + 90 deg 
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Table 19. Results of Parametric Analyses. 


Run 

No. 

Description 

Final 
Position 
Error, Ft 

RMS 
Position 
Error, Ft 

Final 
Altitude 
Error, Ft 

RMS 
Altitude 
Error, Ft 

Ksn 

Nominal Spiral Approach, 3-D 

1 .02 

19.3 

0.34 

1.41 


Initial Crosstrack Error, 1 nm East 

1.02 

1667.7 

0.34 

1.40 


Initial Crosstrack Error, 1 nm West 

1.05 

1681.7 

0.36 

1.46 

1004 

Initial Along-Track Error, 1 nm North 

1.02 

1866.2 

0.34 

1.37 

1005 

Initial Along-Track Error, 1 nm South 

1.03 

2003 . 1 

0.35 

1.45 

1006 

Initial Groundspeed Error, +20 Kt 

1.05 

66.7 

0.36 

1.47 

3007 

Initial Groundspeed Error, -20 Kt 

1.05 

52.1 

0.35 

1.43 

1008 

Initial Altitude Error, +300 Ft 

1.02 

19.3 

0.34 

48.36 

1009 

Initial Altitude Error, -300 Ft 

1.02 

19.3 

0.34 

48.38 

1010 

Initial Track Error, +30 Deg 

1 .03 

139.3 

0.35 

1.40 

1011 

Initial Track Error, -30 Deg 

1.06 

140.1 

0.36 

1.39 

1012 

Wind Speed, +15 Kt 

1.06 

19.4 

0.34 

1.77 

1013 

Wind Speed, -15 Kt 

0.43 

18.8 

0.46 

1.64 

1014 

Wind Direction, +90 Deg 

0.99 

19.1 

0.27 

1 .63 

1015 

Wind Direction, -90 Deg 

12.54 

19.7 

0.71 

1.77 

1016 

Wind Speed Shear, +15 Kt 

0.96 

19.1 

0.30 

1.44 

1017 

Wind Speed Shear, -15 Kt 

0.09 

19.4 

0.52 

1.43 

1018 

Wind Direction Shear, +90 Deg 

1.00 

19.3 

0.33 

1 .42 

1019 

Wind Direction Shear, -90 Deg 

2.33 

19.3 

0.25 

1.43 

1020 

Guidance Parameter, C^ = 0.001 

0.34 

31.9 

0.04 

1.42 

1021 

Guidance Parameter, Cp = 0.010 

0.06 

18.9 

0.04 

1.42 

1022 

Guidance Parameter, Cijr = -1.0 

0.06 

19.8 

0.04 

1.42 

1023 

Guidance Parameter, C^ji = -3.0 

0.16 

20.4 

0.04 

1.42 

1024 

Guidance Parameter, K^ = -0.1 

0.11 

19.7 

0.04 

1.39 

1025 

Guidance Parameter, Ky = -0.3 

0.19 

18.5 

0.07 

1.42 

1026 

Guidance Parameter, Kj^ = 0.5 

0.11 

19.3 

0.06 

1.98 

1027 

Guidance Parameter, Kj^ = 2.0 

0.11 

19.3 

0.03 

1.15 

1028 

Guidance Parameter, 9^, = 45 Deg 

0.12 

19.3 

0.04 

1.42 

1029 

Guidance Parameter, 9^ = 180 Deg 

0.11 

19.3 

0.04 

1.42 

1030 

Guidance Parameter, W* = 30 Kt 

0.34 

19.1 

0.32 

1.15 

1031 

Guidance Parameter, W* = 0 Kt 

2.49 

20.0 

1.73 

1 .95 

1032 

Guidance Parameter, T]* = 180 Deg 

2.06 

20.2 

0.27 

2.35 

1033 

Guidance Parameter, Tj* = 0 Deg 

32.02 

21.1 

5.90 

1 

1.90 
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Figure 41 . Initial Along-track Error of 1 nm. 
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Figure 45. Constant Wind Direction Error of -90 deg. 
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SECTION 9 

SUMMARY AND CONCLUSIONS 


Guidance and navigation requirements of a spiral descent approach concept 
for commercial VTOL operations have been Investigated. The development of the 
spiral descent as a standard IFR approach procedure is an attractive and feasible 
means of satisfying several requirements for safe, efficient and independent VTOL 
operations. Although the spiral descent was first conceived as a means of providing 
independent IFR operations for VTOL aircraft at conventional airports, the technique 
is equally applicable to any heliport with airspace restrictions due to other traffic, 
noise constraints or obstructions. 

The spiral descent retains most of the advantages of a vertical descent, but 
requires less power and fuel, maintains forward airspeed and controllability, and 
avoids the vortex ring state. The problem with the spiral descent operation is that 
it presents a difficult guidance and navigation task, particularly in a wind. The 
steady-state turn introduces an additional integration into the guidance and control 
loop, making it more difficult to stabilize. 

The effect of the wind on the spiral descent is to require a continuous varia- 
tion in the rotorcraft bank angle and heading rate to remain on the desired path. The 
choice of maintaining constant airspeed rather than constant groundspeed during the 
descent was made for several reasons: the descent airspeed can be selected for fuel 
economy; large airspeed variations could place the rotorcraft in unsafe flight regimes; 
and continually-varying longitudinal accelerations would be necessary to maintain 
constant groundspeed. 

The basic concept adopted for the spiral approach guidance is that of a 
feedback system to keep the rotorcraft close to a prespecified nominal path. To 
simplify the analytical formulation, particularly for 3-D guidance, the spiral turn 
angle was used instead of time as the independent variable. For 4-D guidance, the 
nominal time must be computed as a function of the turn angle, which was shown to 
require the solution of an elliptic integral of the second kind. A truncated series 
expansion of the elliptic integral, keeping only the first two terms, was found to 
provide an accuracy of better than 0.3% in approximating the nominal time calculation. 
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The nominal spiral is based on constant' radius, constant airspeed, constant 
rate of descent, and steady predicted winds. The controlled states are radial position, 
heading, airspeed, altitude and time (for 4-D guidance); the guidance outputs are 
commanded values of bank angle, longitudinal acceleration, and vertical speed. 

For passenger safety and comfort, the guidance commands must be limited to reasonable 
levels, 

A linear perturbation analysis was used to select a satisfactory set of guidance 
gains for preliminary evaluation of the guidance concept. For 4-D guidance, time 
errors can be corrected by airspeed adjustments and/or by small changes in the radius. 
However, the latter technique couples the longitudinal and lateral guidance loops, 
and somewhat destabilizes the 3-D lateral response. Consequently, the feedback 
gain term providing bank corrections to timing errors was omitted from further analysis. 
The resulting guidance gains were used as initial values in the interactive nonlinear 
simulation. An iterative evaluation produced an improved set of gains, along with 
additional modifications to the guidance scheme. 

As expected, the guidance system performance was found to be extremely 
sensitive to errors in the predicted wind used for the nominal calculations. Con- 
sequently, a wind estimator will be an essential element of the guidance system, 
both for 3-D and for 4-D guidance. The updated wind estimates are used to revise 
the nominal calculations of bank angle and heading during the descent. However, 
for 4-D approaches, the nominal time calculation uses the original wind prediction 
to maintain proper spacing between aircraft. 

Two approaches were investigated for providing the onboard wind estimates. 

A Kalman filter was formulated using the linearized state equations with the guidance 
loops closed. The north and east wind components were each modeled as exponentially- 
correlated random processes (first-order colored noises). The two wind estimate 
equations must be augmented to the state equations, giving a 7th order error state 
vector. To implement the estimator, a 7 x 7 error covariance matrix must be 
propagated . 
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In lieu of the Kalman filter approach, a simplified wind estimator was 
developed by using the off-nominal values of bank angle and heading angle to indicate 
errors in the tangential and radial wind components, respectively. The simplified 
estimator requires the specification of a characteristic angle, over which the wind 
estimates are averaged. In a descent through a constant wind, initial errors in the 
estimate will be reduced to 1/e over the characteristic turn angle; in a constant shear, 
the wind estimate will lag the actual wind by approximately the chcmcterlstic angle. 
Simulation results have vei ified the performance of the simplified wind estimator, 
using a characteristic angle of 90°. 

A phase plane analysis was used to analyze the behavior of the bank guidance 
loop for large errors in radial position and/or heading. In such situations, the bank 
angle is constrained and the system response is nonlinear. For initial spiral offsets 
greater than 1200-1500 feet, the bank guidance law will command a continuous maxi- 
mum bank angle, resulting in a tight circular turn which never Intercepts the desired 
path. To ensure spiral capture, the bank guidance law was modified to produce a 
constant nominal intercept angle for arbitrarily large radial errors. An intercept 
angle of 60 degrees was found to be satisfactory with the interactive simulation. 

The transitions from the enroute guidance mode to the spiral, and from the 
spiral to hover and touchdown, are critical phases which must be accomplished 
smoothly. For the present study, these transition paths have been defined as straight 
flight segments tangent to the nominal spiral . For consistency with the spiral itself, 
the independent variable has been selected as along-track distance In place of time. 

The nominal entry transition assumes constant altitude and constant deceleration to 
the nominal spiral airspeed, while the exit transition assumes a constant rate of descent. 
The guidance laws for the transition segments are nearly identical to those during the 
spiral, with the addition of a nominal deceleration term in the acceleration loop. 

The bank guidance loop is modified to ensure capture in the same manner as in the 
spiral. Simulation results have verified the performance of the transition segment 
guidance algorithms. 

Two digital simulation programs were developed for analyzing and evaluating 
the spiral guidance concept. An inhouse, interactive simulation, containing a simple 
model of the vehicle dynamics and no navigation errors, was used to verify the guid- 
ance and wind estimator algorithms and to study guidance parameter adjustments. 
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This simulation was especially valuable In providing a clear understanding of the 
entire spiral guidance problem. 

The second simulation is a modification of Program VALT, and was used to 
evaluate the performance of the complete spiral guidance concept for a typical approach 
to New York City's Kennedy International Airport. The VALT simulation includes both 
the enroute guidance and spiral guidance algorithms, a general model of the deter- 
ministic and random wind inputs, stochastic models of the various navigation sensors, 
and a Kalman filter for generating the navigation system outputs. The VALT simu- 
lation was exercised for a number of parametric variations to evaluate the performance 
of the spiral guidance concept. Initial condition errors can produce significant errors 
in the navigation system due to nonlinear effects, as well as timing errors for 4-D 
guidance. Wind errors give essentially the same results as the interactive simulation; 
the simplified estimator determines the winds with a delay equivalent to the charac- 
teristic turn angle. Variations in the guidance parameters verify the analytical results 
obtained with the linearized system model . 

In summary, the investigation has demonstrated that the spiral guidance tech- 
nique is technically feasible for commercial VTOL operations. Wind errors have a 
very important effect on the guidance system performance. Consequently, a wind 
estimator will be required for the spiral approach to adfust the nominal path appro- 
priately. Finally, the spiral guidance algorithms and models are general, and can be 
applied easily to various rotorcraft and landing sites. 
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SECTION 10 


j 


RECOMMENDATIONS FOR ADDITIONAL RESEARCH 

During the course of this study, several areas have been identified In which 
additional research is needed to further evaluate the feasibility of the spiral approach 
guidance concept for commercial VTOL operations. These subject areas are outlined 
below. 

10.1 SPIRAL GUIDANCE FLIGHT TEST PROGRAM 

A full-scale flight evaluation of the spiral guidance technique should 
eventually be conducted as part of the VALT flight test program. In anticipation of 
this, a detailed flight test plan should be prepared to integrate the spiral guidance 
tests with the other planned experiments. This is especially important in the present 
VALT organization since several research groups comprising a variety of technical 
disciplines will be involved in planning, conducting and evaluating the flight test 
program . 

A preliminary outline of the major elements to be considered in the spiral 
guidance flight evaluation is presented in the following subsections. 

10.1.1 PURPOSE OF TEST 

• Verify conclusions of computer simulation and analytic studies by 
actual flights under automatic control 

• Evaluate ability of human pilot to fly manual spirals using flight 
director commands 

• Obtain pilot evaluation of operational procedures 

• Check for possible disorientation due to vertigo 

• Determine any impact of spiral on passenger comfort 

10.1.2 TEST EQUIPMENT AND FACILITIES 

• VALT CH-47 helicopter test bed including measurements of position, 
altitude, groundspeed, vertical speed, airspeed, heading, bank angle 
and software 

• Attitude autopilot 

• Flight director 
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• Timely wind soundings 

• Ground tracking radar 

• Test pilots - hooded subject plus safety pilot 

• Flight test support facilities 

• Post-flight data reduction 

10.1.3 FLIGHT PROCEDURES 

Fly a large matrix of spirals with following variations: 

• Automatic and manual control 

• Various rates of descent from level to 1500 ft/min 

• 3-D and 4-D guidance algorithms 

• Various test pilots 

n Various wind conditions including shear and turbulence if practicable 

• Various initial conditions to evaluate capture logic 

• Various radii of spiral 

• Various speeds 

• With and without wind estimate or variations of guidance gains 

10.1.4 DATA OBTAINED 

• Ground-based measurements of position and velocity versus time 

• Airborne measurements of airspeed, heading, altitude, attitude and 
accelerations 

• Wind estimates based on balloon soundings 

• Guidance system commands 

• Test pilot Cooper ratings 
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10.1.5 DATA ANALYSIS 


• Reduce fllghl- data with plots and statistics of results 

• Attempt to demonstrate that flight results were predicted by preliminary 
analysis 

• Explain differences between analysis and flight results 

• Compare manual and automatic results 

• Assess practicability of spiral descent as operational technique 

• Recommend procedures, parameter values, limitations, etc., based on 
combined results of analytic and flight evaluations. 

10.2 GUIDANCE FOR SPIRAL EXIT AND LANDING 

Work on spiral descent guidance under the present contract has shown the 
advantages and feasibility of this VTOL approach technique. The spiral enables a 
VTOL aircraft to avoid obstacles and noise sensitive areas, to operate independent of 
CTOL traffic patterns, and continue to maintain reasonable airspeed for increased 
safety. It allows the VTOL vehicle to descend into confined air space with effective 
glide slope angles that are relatively shallow. Exit from the spiral occurs nominally 
at 60 knots at an altitude from 200 to 400 feet and is followed by deceleration and 
landing. The exit and landing phase needs additional study to determine the placement 
of the spiral relative to the landing pad. 

Present studies have been limited to landing along a tangent to the spiral (see 
Figure 49). However, alternate concepts, such as landing inside of the spiral , along 
the spiral, or outside the spiral are feasible and could be advantageous for some 



Figure 49 . Possible Heliport Locations Relative to the Descent Spiral, 
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heliports. A related question to be considered is whether there should be only one 
spiral for several landing pads or a dedicated spiral for each pad. Factors which 
influence these decisions are: 

• Geography of the landing area 

• Available navigation coverage 

• Operational flexibility 

• Magnitude of the necessary maneuvers 

• Possibility of pilot vertigo 

• Difficulty of the guidance and control tasks 

• Location of CTOL air traffic patterns 

• Location of noise sensitive areas 

Guidance of the VTOL vehicle must undergo a transition from an airspeed 
control mode during approach to a hover mode prior to landing. The airspeed control 
mode is characterized by coordinated turns using roll to control the direction of the 
airspeed vector; whereas, the hover mode is characterized by heading control through 
the yaw channel with pitch and roil used to maintain ground position and velocity. 

The transition will normally occur after exit from the spiral during the deceleration 
for landing. To accomplish this transition, a description of the nominal path to each 
available landing pad must be formulated, along with a guidance scheme to maintain 
that path during the transition. The scheme then should be evaluated to determine 
its performance in the presence of winds, guidance errors, and parameters in the 
guidance algorithms. 

10.3 WIND ESTIMATION IN THE VALT FLIGHT PROGRAM 

The present onboard computer software configuration for the Langley VALT 
flight program guidance logic uses a flight-path-following system that seeks to null 
cross-track position error and cross-track ground speed error. The position is measured 
by a ground tracking radar, and the ground speed is derived from the position measure- 
ments. In an operational system, ground speed information would not generally be 
ovailable from the ground. Although onboard inertial navigation systems or Doppler 
radars could be used to provide the ground-referenced velocity, both of these equip- 
ments are relatively expensive compared to more conventional air data and vehicle 
attitude measurements. 
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Flight-path following without ground speed information substitutes heading 
error for cross-track ground speed error. A wind estimate or integration in the feed- 
back path of the cross-track position error is required to prevent a steady-state posi- 
tion offset. Wind estimation is preferable since heading changes introduce a variable 
cross-track wind component, even when the wind velocity is constant relative to the 
ground. The present VALT flight computer software can be modified to utilize the 
same wind estimator developed here for spiral descent guidance. The nominal flight 
path developed for the flight program consists of patched elliptical arcs which can 
degenerate into circles and straight lines. Although the AS! wind estimator was 
developed here for use along circular segments, the concept can be extended to 
handle elliptical arcs. This would give the VALT flight guidance system an additional 
capability for wind estimation without ground speed measurement, and consequently 
would provide more accurate tracking when only air data and attitude information are 
available. This, in turn, would lead to a more realistic assessment of the requirements 
for inertial navigation or Doppler radar. 

Consequently, an evaluation should be made of the performance of the 
present onboard computer guidance logic for the VALT flight program in the 
absence of ground speed information. A wind estimator compatible with the guidance 
logic utilized on elliptical flight paths should be designed using measurements of 
position, attitude and airspeed but without ground speed information. Existing com- 
puter programs could be modified to simulate the wind estimation scheme when operating 
along elliptical arcs. The simulation results would predict the performance comparison 
between the existing VALT flight guidance logic with ground speed information availa- 
ble, as opposed to the modified logic without ground speed information but with the 
new wind estimator. Flight experiments could test the validity of conclusions drawn 
from the computer simulation. 
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